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A material whose electrons are correlated can affect electron dynamics across the interface with 
another material. Such a "proximity effect" can have several manifestations, from order parameter 
leakage to generated effective interactions. The resulting combination of induced electron correla- 
tions and their intrinsic dynamics at the surface of the affected material can give rise to qualitatively 
new quantum states. For example, the leaking of a superconducting order parameter into certain 
Rashba spin-orbit-coupled materials has been recently identified as a path to creating "topological 
superconductors" that can host Majorana particles of use in quantum computing. Here we ana- 
lyze the other aspects of the superconducting proximity effect. The proximity-induced interactions 
are a promising path to incompressible quantum liquids with non-Abelian fractional quasiparticles 
in topological insulator quantum wells, which could also find applications in topological quantum 
computing. We discuss the operational and design principles of a heterostructure device that could 
realize such states. We apply field-theoretical methods to characterize the properties of induced 
interactions via the electron-phonon coupling and Cooper pair tunneling across the interface. We 
argue that bound-state Cooper pairs can be stabilized by the interaction proximity effect inside a 
topological insulator quantum well at experimentally observable energy scales. The condensation of 
spinful triplet pairs, enabled by the Rashba spin-orbit coupling and tunable by gate voltage, would 
lead to novel superconducting states and fractional topological insulators. 



I. INTRODUCTION 

Interfaces between materials with radically different 
electronic ground states are increasingly promising plat- 
forms for engineering novel quantum states. The quest 
for elusive Majorana particles has been driving the latest 
wave of interest in such interfaces. Majorana quasiparti- 
cles were first predicted to exist as special Andreev bound 
states at the boundaries or inside vortex cores of spinless 
weakly-coupled p x + ip y superconductors^, and similarly 
at the terminal points of quantum wires^. However, it 
seems more practical to seek analogous Majorana par- 
ticles in superconducting states with time-reversal (TR) 
symmetry that can be artificially created inside mate- 
rials with a strong spin-orbit coupling. Such materials 
are not natural superconductors, but pairing can be in- 
duced in them via the "proximity effect" , by placing them 
in contact with a conventional superconductor 5 . Nu- 
merous theoretical proposals of proximity effect devices 
utilizing surfaces of the bulk topological insulators^— 
(TI), films^ and quantum wires£~— have been made and 
accompanied by microscopic model calculations^^— . 
The experimental search for Majorana particles based 
on those ideas 2 ^— and the exploration of proximity 
effects 2 ^— are also gaining momentum. 

Consider an interface between a correlated material 
C whose ground state spontaneously breaks some sym- 
metry (such as a superconductor) and an uncorrelated 
material UC whose ground state is naturally disordered 
(such as a metal or band- insulator). The term "prox- 
imity effect" is usually used to describe the leakage of 
the order parameter across the interface from C to UC. 



The material C presents an explicit symmetry-breaking 
perturbation to the electron dynamics inside UC, so that 
some electron correlations of the kind found in C will in 
most cases appear near the interface in UC (if its sus- 
ceptibility to the perturbation is finite). In this manner 
correlations can be imposed on electrons in UC where 
they do not naturally occur. This phenomenon is be- 
hind all proposals to obtain Majorana particles in TIs 
and similar systems. 

In this paper we qualitatively analyze other ways in 
which a correlated material C can induce correlations 
among electrons in a naturally uncorrelated material UC. 
Our interest is the "interaction proximity effect" , or the 
renormalization of electron interactions in UC by their 
coupling to the degrees of freedom in C. We will specifi- 
cally consider an interface with a conventional supercon- 
ductor (SC) and explore two sources of effective inter- 
actions in the two-dimensional electron system (2DES) 
formed near the interface in the material UC: the direct 
coupling of the 2DES electrons to the SC's phonons, and 
Cooper pair tunneling from the SC to the 2DES. The 
second mechanism incorporates the usual "proximity ef- 
fect" . We will describe the proximity-induced interac- 
tions in 2DES using interaction potentials and effective 
actions, and discuss their ability to generate new states of 
matter with spontaneously broken symmetries or topo- 
logical order. 

The current motivation to study the proximity-induced 
interactions comes from the possibility of obtaining 
fractional topological insulator states in TI quantum 
wells 38 ' 39 . Consider a heterostructure in the FigJT] which 
consists of a TI quantum well sandwiched between a bulk 
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FIG. 1: (color online) The heterostructure device that embeds 
a topological insulator quantum well (TI) between a conven- 
tional superconductor (SC) and a conventional insulator (I). 
The chemical potential fj, in the quantum well can be con- 
trolled by applying a gate voltage V g , and more efficiently so 
if a tunneling barrier (B) is inserted at the SC-TI interface. 
A schematic band-structure of this device is also shown. Bulk 
band edges are drawn by thin lines, while the two-dimensional 
band edges of the quantum well are represented by red lines. 
These two-dimensional bands are created by the hybridiza- 
tion of surface states that would be gapless and topologically 
protected in a bulk three-dimensional TI. Their surface hy- 
bridization bandgap is smaller than the bandgap of a bulk TI. 
The blue dotted lines are the two-dimensional sub-gap Cooper 
pair bands created by the interaction proximity effect. 

superconductor (SC) and a conventional insulator (I). 
The TI-SC interface can be bridged by a tunneling bar- 
rier (B), but this may not be necessary. This setup allows 
one to control the chemical potential inside the insulating 
TI quantum well by the gate voltage V g and hence drive 
quantum phase transitions at low temperatures. We will 
consider quantum wells made from "strong" TI materials 
such as Bi2Se3 or Bi2Te3. Electrons inside the quantum 
well have the spin a z and a two-state "orbital" degree of 
freedom r z . If we neglect their coupling to the rest of 
the heterostructure, then the basic aspects of their two- 
dimensional dynamics are captured by the single-particle 
Hamiltonian written in terms of the spin a x,y,z and or- 
bital t x ' v ' z Pauli matrices: 

H = V -i(a x p)r z + At x -h , (1) 

where v is the strength of the Rashba spin-orbit coupling, 
and A is a bandgap formed by the hybridization between 
the helical Dirac quasiparticles from the opposite sur- 
faces t z — ±1 of the TI quantum well. The proximity- 
induced interactions among these electrons of any spin 
can create two-dimensional Cooper pairs in several chan- 
nels. Particularly interesting are the inter-orbital spin- 
triplet pairs with a finite spin projection. Such pairs 
are turned into low energy excitations by the Rashba 
spin-orbit coupling in the TI. They can live at ener- 
gies inside the fermionic bandgap and "condense" by 



the application of a gate voltage. This "condensate" is 
expected to host a TR-invariant vortex lattice, shaped 
by the strong SU(2) "magnetic" flux of the spin-orbit 
coupling 3 ^ (strictly speaking, only the translation sym- 
metry is spontaneously broken in this state, and only 
at zero temperature 4 ^ - — ) . The phase transition to a 
non-superconducting state of triplets is first-order at suf- 
ficiently low temperatures and originates in the zero- 
temperature vortex lattice melting by quantum fluctua- 
tions. The ensuing insulating state is strongly correlated 
and a candidate for novel fractional incompressible liq- 
uids with non-Abelian quasiparticle a 43 ' 44 . All these cor- 
related states of triplet Cooper pairs can coexist in the 
TI with a singlet condensate. 

When a gate voltage V g is applied across the device 
to attract electrons from the SC to the TI, the band- 
structure of the bulk spectrum assumes the qualitative 
shape shown in the FigJTJ If the TI material were thick, 
its surfaces would host a topologically protected band of 
massless Dirac quasiparticles with a helical locking be- 
tween their spin and momentum. However, the quan- 
tum well is thin enough to hybridize these surface states 
on its opposite sides, so a two-dimensional bandgap 2A 
opens and the surface Dirac quasiparticles become mas- 
sive. The hybridized surface states still live at energies 
well below the TPs bulk bands, and their bandgap 4 ^ can 
be of the order of A - 10 - 100 meV. The SC's spec- 
trum consists of Bogoliubov quasiparticles whose (pair- 
ing) gap Ao ~ 1 meV is set by the critical temperature 
T c ~ 10 K. Charge is very mobile inside the SC, so this 
layer (the material C) provides a particle reservoir for the 
TI quantum well (the material UC) and fixes the chem- 
ical potential fi shown by the dashed line in the Figffl 
The same chemical potential must be established in the 
TI because the interface and the tunneling barrier do not 
prevent the transfer of particles. The gate voltage then 
controls the relative energy between the TPs and SC's 
electrons, or equivalently the position of the TPs chemi- 
cal potential relative to its bands. Since only a fraction 
AV~ ~ V g (lB + d)/(lB+d+h) of the applied voltage affects 
the quantum well, it is desirable that the insulator thick- 
ness l\ be as small as possible without jeopardizing the 
insulating behavior of I. Note that the two-dimensional 
character of the TIs bands is protected in the range of 
energies n ± Ao where there are no bulk SC states that 
could resonantly couple to the TI across the interface. 

This device has an interesting application described 
earlier only if the chemical potential can be effectively 
moved in the vicinity of a surface band edge in the TI. A 
generic solution would be to engineer a sufficiently thick 
tunneling barrier that can support voltage drops of the 
order of the bulk bandgaps (~ 1 eV) without suffering 
an electrical breakdown. This would, unfortunately, ruin 
the proximity effect. Instead, a barrier Zb ~ 10 A thin 
could be theoretically used to tune the chemical poten- 
tial in the TI within a range of a few or at most ten 
meV (using AlGaAs as an example) . This would be quite 
enough to drive the subtle quantum phase transitions of 
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triplet Cooper pairs, but still requires that the materi- 
als' interfaces in the heterostructure naturally put the 
TI's chemical potential within about 10 meV from one 
of its hybridized band edges. Therefore, a careful choice 
of materials, doping, etc. may be needed to obtain the 
adequate work functions and other interface properties. 
The interaction proximity effect might still be felt a short 
distance away, helped by the fact that an insulating 
barrier cannot screen the charge of the SC's electrons 
and atoms. A tunneling barrier could be even unneces- 
sary, and removing it would be the best for proximity 
effect. Namely, the insulating TI can itself play the role 
of a tunneling barrier for the purposes of creating inter- 
orbital spin triplet Cooper pairs. Its desired thickness of 
the order of d ~ 10 A is likely sufficient to tune the chem- 
ical potential across the low-energy landscape of possible 
correlated phases. This characteristic range of energies 
for correlated states is bounded by the SC's critical tem- 
perature, which is at best of the order of 1 meV. 

Here, however, our goal is not to investigate the mi- 
croscopic device design, but rather the device's funda- 
mental physical feasibility and principles of its opera- 
tion. We will particularly show that a sufficiently strong 
interaction proximity effect can create a triplet bound- 
state Cooper pair in the TI, whose band lives below the 
fcrmionic quasiparticle continuum. The bound-state is 
made possible by the two-dimensional geometry of the 
quantum well^ and its insulating ground state. Further- 
more, the triplet bound states with large net momen- 
tum can gain energy from their Rashba spin-orbit cou- 
pling and become even more advantageous than singlets. 
Putting the chemical potential in their low-energy band 
creates a triplet "condensate" (in the mean-field approxi- 
mation), which in this case hosts the TR- invariant vortex 
lattice mentioned earlier. 

The physics portrayed here should be contrasted with 
other examples of proximity effects and two-dimensional 
superconductivity. First of all, the proximity effect ap- 
plied to a 2DES with no orbital degree of freedom, such 
as the surface of a bulk TI, can hardly produce a triplet 
Cooper pair of sufficiently small size whose spin-orbit 
coupling could turn it into a competitor to singlets for 
condensation. Furthermore, no Cooper pair could exist 
as a bound state when the fermion quasiparticle spec- 
trum is gapless. In contrast, the system in the Fig.l 
provides an orbital degree of freedom (surface index) to 
electrons, confines them in the z direction to allow small 
inter-orbital triplet pairs, and maintains an energy gap in 
the quasiparticle spectrum which can be partially filled 
by bosonic excitations. This is the basis for the exis- 
tance of novel correlated states in this system. Should 
a non-Abelian fractional TI be realized in this system, 
it might provide a more robust platform for topological 
quantum computation than the Majorana quasiparticle 
route. Namely, Majorana quasiparticles remain topolog- 
ically protected only while being kept a large distance 
apart, while the non-Abelian quasiparticles of a fractional 
TI are naturally gapped and keep their mutual fractional 



statistics at finite distances. 

A Cooper pair bound state is not normally encountered 
in condensed matter systems. Its existence here is tied 
to the "band-insulating" state of the two-dimensional 
TI, where any amount of (induced) attractive poten- 
tial between particles creates a bound-state^—. This 
pre-formed Cooper pair still has quite a large size, so it 
quickly loses its identity in the correlated finite-density 
ground states. The main consequence of the pair binding 
energy is that strongly correlated quantum states at low- 
densities are possible, which is particularly important 
for the goal of obtaining fractional incompressible liquids 
(where the density of particles should be of the order or 
smaller than the density of flux quanta). In other words, 
the insulator-superconductor phase transition (of spin- 
ful triplets) is not merely a pairing transition, but a less 
conventional one: first-order, or belonging to the bosonic 
mean-field universality class if second-orde r 38 i 50 ~ — . One 
of our main results is that the Cooper pair binding en- 
ergy, which sets the scale for the correlated phases of 
interest, need not be much smaller than the SC's critical 
temperature scale T c ~ 10 K. Thus, the novel interesting 
correlated states should be experimentally accessible in 
the system of this kind. 

Even though the interaction proximity effect need not 
result with novel correlated states of matter in arbitrary 
interfaces between two materials, we will study it quite 
generally. Our analysis will not make any concrete as- 
sumptions about the degrees of freedom available in the 
affected 2DES, or their dynamics. We will only determine 
the ways in which the SC's degrees of freedom renormal- 
ize or induce the interactions among the 2DES electrons. 
A conventional superconductor has two important de- 
grees of freedom: its lattice phonon modes, and BCS- 
paired electrons. 



A. The superconducting proximity effect 

Our discussion will focus on the interfaces between con- 
ventional superconductors and TIs, following the current 
experimental efforts^— . There are two sources of ef- 
fective attractive interactions among electrons in the TI 
near the interface: the phonon mechanism and virtual 
Cooper pair tunneling. The former is the same micro- 
scopic BCS mechanism that is responsible for pairing in 
the SC material. Any charge-carrying excitation of the 
TI's surface electrons can displace the atoms of the SC 
crystal near the interface, and thus attractively interact 
with another TI's excitation via emission and absorption 
of the SC's phonons. The Coulomb origin of this inter- 
action makes it spin-independent and capable of creating 
triplet Cooper pairs in principle. Its range is either short 
due to screening, or blunted by geometry (we describe 
this in detail in the section III D[) . This attractive inter- 
action is further aided by the phonons of the remaining 
TI's environment. 

The second source of effective interactions is the dy- 
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namics of the SC's Cooper pairs. Electrons from the 
SC can gain kinetic energy by tunneling into the TI and 
back, provided that they can maintain their pairing cor- 
relation while being inside the TI. Therefore, tunneling 
can dynamically generate pairing forces in the TI, which 
are short-ranged due to the SC's Meissner effect. These 
dynamic forces reflect the SC's pairing and thus are spin- 
dependent, but we will show that they nevertheless also 
contribute to the possibility of triplet pairing. Cooper 
pair tunneling also produces the conventional proximity 
effect, a direct imprint of the SC's s-wave order parame- 
ter across the interface. This is important in a metallic 
surface 2DES, but becomes only a small perturbative ef- 
fect in the TI quantum wells of our interest (the explicitly 
induced pairing gap is ~ Aq/A « A ). 

The total generated attractive interaction in the TI 
has to be sufficiently strong to overcome the Coulomb 
repulsion and form Cooper pairs. Fortunately, it seems 
that the needed conditions for pairing and superconduc- 
tivity can be achieved in the gated device from the FigfT] 
Pairing in the inter-orbital spin-triplet channels, between 
two electrons of equal spin on the opposite quantum 
well surfaces, is also feasible. Namely, the Rashba spin- 
orbit coupling lowers the energy of triplet pairs, and the 
Coulomb forces that the phonon mechanism is based on 
are unscreened or poorly screened across the quantum 
well. Some in-plane screening inside the TI is always 
provided by the SC's electrons via the image-charge ef- 
fect, but screening across a thin quantum well cannot be 
effective even in the proximity-induced superconducting 
state because the TPs Cooper pair condensate is funda- 
mentally two-dimensional so its wavefunction is rigidly 
determined in the z-direction by the quantum well con- 
finement. The dynamics of some TPs electron's remains 
two-dimensional so they can form bound-state Cooper 
pairs. 



B. The contents and main results 

The first part of the paper, section [Til analyzes the 
phonon-mediated interactions. We will first review in 
the section Hi Al the BCS theory prediction that the pair- 
ing gap is not a function of the underlying Fermi surface 
size in two-dimensions. This is in stark contrast to three- 
dimensional superconductors. We will show in the sec- 
tion [TIB] that this result naively extrapolates all the way 
to the zero density, where the Fermi surface is shrunk and 
disposed of in favor of a band-insulating ground state. 
We will recall another well-known result that a bound- 
state always exists in two-dimensional attractive poten- 
tials, and show that the binding energy has the same 
dependence on the electron-phonon coupling as the 2D 
superconductor's pairing gap. This calculation is sim- 
ple, but ignores the phonon dynamics. We will redo the 
calculation in the section III CI using the Bethe-Salpeter 
equation to fully take into account the phonon dynam- 
ics, and argue that the Cooper pair bound-state remains 



stable despite retardation effects. Finally, we will an- 
alyze the fate of the bound-state against the Coulomb 
repulsion in the section III Dl and show that it survives 
when the electron-phonon coupling is large enough. The 
bound-state criterion is a slight geometric modification 
of the analogous criterion for the existance of supercon- 
ductivity in the 3D bulk SC material. Based on this we 
argue that a sufficiently thin quantum well (and tunnel- 
ing barrier) should allow the phonon-mediated pairing 
in realistic SC-2DES interfaces. The binding energy esti- 
mates in this paper do not include any effects of the spin- 
orbit couplings. Such effects have been taken into ac- 
count elsewhere^, and only help the triplet Cooper pairs 
to form a two-dimensional band within the fermionic 
bandgap. 

In the second part of the paper, section IIIII we con- 
sider the tunneling effect of Cooper pairs on the 2DES. 
We model the interface simply as a tunneling barrier and 
derive the effective action of the 2DES electrons by in- 
tegrating out the SC's electron degrees of freedom using 
a path-integral. The section HlI Al explores the simplest 
consequences of tunneling where only the shortest-range 
interactions are induced. The effective action contains 
a direct symmetry-breaking coupling between the 2DES 
electrons and the SC's order parameter (the conventional 
proximity effect), and also an effective interaction term. 
The latter contributes to the same type of correlation as 
the conventional proximity effect. In the following section 
IIIIBI we analyze the tunneling-generated interactions in 
any available inter-orbital channels, which also includes 
spin-triplet channels. We find that such attractive inter- 
actions are generally induced by quantum fluctuations, 
but probably too weak to bring about any new physics on 
their own. They may, however, help the phonon mech- 
anism to generate triplet pairing. Aspects of Coulomb 
interactions between the 2DES and SC electrons are cov- 
ered in the section Hi Dl 



II. PHONON MECHANISM 

Here we focus on the phonon-mediated proximity effect 
in a band-insulator quantum well that touches a conven- 
tional superconductor. We will argue that the phonon 
mechanism can give rise to bound-state Cooper pairs in 
the two-dimensional band-insulator, whose binding en- 
ergy scale may be large enough to observe. This scale 
determines the critical chemical potential (away from 
the electron's conduction or valence band) for the zero- 
temperature superconductor-insulator transition, as well 
as the typical critical temperature deeper in the super- 
conducting phase. The analysis has four parts that grad- 
ually build a full picture of the proximity phonon mech- 
anism and eventually take the Coulomb repulsion into 
account. 

The coupling between electrons and phonons was stud- 
ied in great detail for a long time^. The minimal model 
of electron and phonon dynamics, which contains all im- 



5 



portant ingredients for conventional superconductivity, is 
given by the Hamiltonian 



H = d a r 



1 

2m 



(V^)(V^)-/i^ 



+ug^Jj a Vcf) + 1 f c(> 2 + u 2 ( V0 V 1 



(2) 




FIG. 2: Phonon-mediated scattering between two electrons. 



in the continuum limit. The real vector field <f> describes 
the "displacement density" of atoms in the crystalline 
lattice, and <f> is its canonical conjugate. ip a is the elec- 
tron field (annihilation) operator. This Hamiltonian is 
valid at energy and momentum scales where the disper- 
sion of sound is approximately determined by a single 
parameter, the sound velocity u ~ 10 3 m/s. The point- 
group symmetries forbid the minimal linear electron den- 
sity coupling to transverse phonons. Any transverse cou- 
pling must involve either scalar combinations of <p or elec- 
tron currents, so it can be neglected perturbatively due 
to the smallness of its coupling constant or electromag- 
netic relativistic effects. The electron G and longitudi- 
nal phonon D propagators are given by the time-ordered 
ground-state expectation values: 

G ffff , = -i(0|TVv(z M )^,«)|0) (3) 

D(x^x^) = -z U 2 (O|T[V0(^)] [V0t(^)] | ) , 

and their Fourier transforms are: 



G a a' {Pfj 



lu - E p + iO+sign(E p ) 

u 2 \ P \ 2 
u 2 - u 2 \p\ 2 + i0+ 



(4) 



The clcctron-phonon coupling g has been associated with 
the vertex. We denote by = (r, t) and p^ = (p, uj) 
the space-time positions and momenta respectively (us- 
ing the units h = 1), and E p = |p| 2 /2m— /i is the electron 
energy relative to the chemical potential y,. 



A. Two-dimensional BCS superconductors 

The BCS theory of superconductivity in metals relies 
on several approximations to predict the properties of the 
superconducting state. First, the electron-phonon cou- 
pling g is assumed to be small, so vertex corrections can 
be neglected. Second, phonon-mediated energy transfers 
between interacting electrons can be neglected on aver- 
age. This is justified in metals with sufficiently large 
electron density: typical momentum transfers among the 
dynamically active electrons near the Fermi surface are 
of the order of Fermi momentum pf = \/2mE{ = mvf, 
where V{ ~ 10 6 m/s 3> u, while the maximum energy 
transfer by phonons is of the order of the Debye frequency 
ojd ~ upf <C E{. The electron-electron scattering process 
represented by the Feynman diagram in FigJ5] can then 



be approximated by an attractive interaction via a static 
instantaneous contact potential of strength Ao = g ■ 

U(p,oj)=g 2 D(p,u J )=g 2 IP ' » -A , (5) 

LU Z — U z \p\ z 

because \u\ <C u|p| ~ wd for typical phonons. In the 
last mean-field-like approximation, the pairing gap is de- 
termined by assuming that all electrons near the Fermi 
surface contribute equally to superconductivity as long 
as they are capable of exchanging energy up to ojd with 
other electrons. The pairing gap Ao at zero temper- 
ature is accordingly calculated from the self-consistent 
condition^: 



Ao f d d p Hu< L) 



\E n 



(2n)° 



= 1 



E l 



A 2 



(6) 



and leads to the well-known result in three-dimensional 
(d = 3) metals: 



2wd exp 



2tt 2 



X mp f 



< WD 



(7) 



The same energy scale also determines the critical tem- 
perature T c ~ A . 

Consider now a hypothetical two-dimensional weak- 
coupling superconductor. The pairing gap is now deter- 
mined from ^) with a two-dimensional electron-phonon 
coupling A: 



A f d 2 p 9(uj b ~\E p \) A 27tp f 



2 J ( 2 ^) 2 Je 2 +A 2 d 2 (2n) 2 v f J ^ 



de 



X 2D 



Therefore, 



i 2tt 
A 2D ~ 2w D exp -- — 
Am 



(9) 



in two dimensions. The striking difference from ([7]) is 
that there is no dependence on the Fermi momentum pf. 
In three dimensions it is important to have a large Fermi 
surface with a large pf, otherwise Ao and the correspond- 
ing critical temperature can be exponentially suppressed 
to an unobservable small value. In contrast, the pairing 
gap and critical temperature in two dimensions do not 
depend on the electron density as long as the approxi- 
mations behind this simple model are valid. Instead of 
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extrapolating this result down to zero density, we will 
consider more appropriate models in the following sec- 
tions. 

The main system of our interest is a quantum well in- 
terfaced with a conventional superconductor. The well 
could resemble a BCS-paired two-dimensional metal if a 
fairly large gate voltage were applied to create a dense 
two-dimensional electron gas in it. A very thin quantum 
well is not likely to contribute significant longitudinal 
phonons of its own, and we will not consider its trans- 
verse phonons either. The latter could be minimally but 
weakly coupled because the inversion symmetry is ab- 
sent in the z-direction. The superconductor's phonons 
may anyway be the most likely ones to have a significant 
coupling to the electrons in the quantum well. This cou- 
pling A is surely smaller, but need not be much smaller 
that the superconductor's coupling adjusted for the dif- 
ferent dimensionality of the quantum well. We will there- 
fore loosely estimate A ~ Ao<i _1 , where d is the quantum 
well thickness, or A ~ \o(d + Zb) _1 if a tunneling bar- 
rier is present as in the FigfT] This is expected to be 
valid when the quantum well and tunneling barrier are 
very thin. The critical temperature of the ensuing two- 
dimensional superconductivity can be of the same order 
as in the three-dimensional superconductor if the tunnel- 
ing barrier is removed, but even a more realistic reduction 
of T c by an order of magnitude or two can be experimen- 
tally observable. 



B. A simple model of Cooper pair bound states 

The previous section has suggested that a supercon- 
ducting state can survive in two-dimensional metals at 
any density of electrons. What if the density is lowered to 
zero? Is it possible to have a phase transition between a 
band-insulator and a paired superconductor across which 
the fermion excitation gap remains finite? 

The renormalization group argument yields a pos- 
itive answer to this question in two dimensions, for 
any strength of attractive interactions^. This is a 
consequence of the fact that any short-range attrac- 
tive potential always produces a two-body bound state 
in two dimensions 4 ^. A two-dimensional Cooper pair 
"molecule" may be very large and fragile, but it always 
has a finite binding energy. If their Bose condensation is 
arranged from the band-insulating state, the second order 
phase transition belongs to the mean-field bosonic uni- 
versality class and hence differs from the standard pairing 
BCS transition. The normal state adjacent to this tran- 
sition is a "pseudogap" bosonic Mott insulator, since its 
lowest energy excitations are gapped Cooper pairs ly- 
ing inside the fermion bandgap. In the most mundane 
circumstances, the Mott insulator is separated from the 
band-insulator in the phase diagram only by a crossover. 
However, if the transition involves triplet Cooper pairs, 
or inter- valley pairing, interesting correlated Cooper pair 
insulators that break some symmetry become a possibil- 



ity in two dimensions. In the case of a vortex lattice 
state predicted to result from triplet superconductivity 
in topological insulator quantum wells, the above sec- 
ond order transition is preempted by the first order vor- 
tex lattice melting that potentially stabilizes a fractional 
topological insulato r 43 i 44 . 

It is this range of interesting possibilities that moti- 
vates us to explore the bound-state Cooper pairs. The 
binding energy sets the temperature scale below which 
the correlated insulators and low-density superconduc- 
tors can be observed in the quantum well. Our main goal 
is to estimate this energy scale and determine if it can 
be reasonably large. For simplicity, we will consider only 
the pairing of standard non-relativistic electrons. More 
interesting pairing in the presence of a strong Rashba 
spin-orbit coupling also reveals a finite binding energy^. 

We will here without any justification neglect energy 
transfers by phonons in two-dimensional band-insulators, 
and approximate the phonon-mcdiatcd scattering be- 
tween electrons by a short-range potential with strength 
A. After getting a crude estimate of the pair binding en- 
ergy, we will consider the effect of energy transfers in the 
next section. The effective second-quantized Hamilto- 
nian in the quantum well derived from the phonon mech- 
anism without energy transfers is: 



H = 



d r 



( 10 ) 

We are suppressing any orbital indices that the electron 
fields may have, and allowing a = a' if an orbital degree 
of freedom exists. The chemical potential /i has been 
replaced by an excitation gap eo > appropriate for a 
band-insulator. Since there are no particle excitations in 
the ground-state, we can solve the two-body problem us- 
ing the equivalent single-quantized Schrodinger equation: 



2m 2m 



r 2 ) *(ri,r 2 ) = E*(n,r 2 ) . 

(11) 

By separating the center-of-mass from the relative mo- 
tion *(n, r 2 ) = if) cm {R)e tW if)(r), where R = (n + r 2 )/2 
and r = ri — r 2 = (r, 9) in cylindrical coordinates, we 
obtain the reduced radial Schrodinger equation in the 
center-of-mass frame for ip(r); 

11 d ( di\)\ I 2 , J(a-r) , , . 

5-( r lT +~^- A 2 V> = #- (12) 

mr or \ or J mr z 7ror 

The reduced mass of the relative motion is taken into ac- 
count, and e = E — 2eo is the Cooper pair binding energy 
if negative. At this point we have introduced an ultra- 
violet cutoff length a and regularized the Dirac delta- 
function in the contact potential: 



S(t) 



lim 



d{a - r) 



(13) 



This is necessary in order to solve the equation (|12p . The 
scale a is defined by the momentum cutoff A = a -1 for 
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the phonon propagator in (j3]). The acoustic pho non- 
mediated interaction between electrons cannot be sig- 
nificant at electron separations smaller that the lattice 
constant, so we expect a to be of the order of lattice 
spacing. 

This is a simple and standard problem in quantum 
mechanics. The solution for r > a is expressed in terms 
of the Bessel functions: 



so that: 



(Vr > a) r))(r) = C 1 Ji{k'r) + C 2 Yi{k'r) 



k 1 = 



The energy e of a bound state is negative, so that the 
parameter k' = in is imaginary and the Bessel functions 
exhibit an exponential dependence on r. Their asymp- 
totic behavior in the r — > oo limit 



Ji(inr) 
Yi (inr) 



1 



2nnr 



->• I 



1 



2~KKT 



e^ l e Kr 



implies that we must choose Ci — iC\ in order to have 
a normalizable wavefunction. Inside the area affected by 
attraction the wavefunction is oscillatory: 



(Vr<a) ^(r)=AJ t {kr) , k = _ | e | 

We must match the wavefunction values and derivatives 
on the left r = a — + and right r = a + + side of the 
potential step. We will do that only for I = because 
the existance of bound states is not guarantied when the 
angular momentum I is finite and a is small enough (def- 
initely the case if mX < it). We shall further anticipate 
that the binding energy |e| satisfies 



/ta -C 1 



k = vmfe] 



(16) 



because a is small and the attractive potential A is weak. 
The approximate form of the wavefunction (|14l) in this 
limit is: 



ip(r) > — log (nr) 

7T 



(17) 



so that the boundary conditions at r — a for I — are: 



iC 

AJ (ka) — — logfwa) 
^(j_ 1 (ka)-J 1 (ka))= — 



(18) 



Dividing these two equations yields an energy quantiza- 
tion condition: 



_ Jo(fca) _ 
J\(ka) 

By the approximation Q16p. 



ka log(«;a) . 



ka = 



mX mX 
(nay w \ 

7T V 7T 



(19) 



(20) 



K 

m 



ma' 



■ exp 



7T J Q (y/m\/-K) 



mX J^y/mX/ir) 



mA«l 1 

> ? exp 



47T 

mX 



(21) 



The final result applies in the limit na <C mX <C 1 and 
should be compared with the pairing gap © of the two- 
dimensional BCS superconductor. The dependence on 
the electron-phonon coupling A is essentially the same, 
which justifies the earlier naive assumption that ^ can 
be extrapolated to very low densities. The difference be- 
tween the exponents of © and (f2"Tj) can be attributed 
to the effective mass reduction in relative motion, which 
is not taken into account in ([§]). The Cooper pair size 
£ ~ k^ 1 a is relatively large. 

The energy scale l/(ma 2 ) obtained here can be much 
larger than wd , but depends on the microscopic details at 
the cutoff scale which we cannot model accurately. This 
result has only a qualitative significance. It reflects the 
fact that pairing is much more efficient in two dimensions 
than in three dimensions. The strict weak-coupling limit 
of pairing does not even exist in two-dimensional sys- 
tems at very low densitie a 51 ' 52 , and superconductors with 
(quasi) two-dimensional dynamics tend to have higher 
critical temperatures (MgB2 is an example). 

An even more important point to emphasize is that 
our two-dimensional system of interest is not isolated as 
the above calculation indirectly assumes, but coupled to 
a bulk superconductor whose spectrum has a continuum 
of extended states at energies more than Ao away from 
the chemical potential /i. Only within the energy range 
2Ao can the dynamics in the quantum well be regarded 
as two-dimensional. We will see in the following section 
that this cuts off the Cooper pair binding energy down 
to no more than the bulk superconductor's pairing gap 
scale A : 



E exp 



An 
mX 



En ~ An 



(22) 



Nevertheless, we will argue in the section III Dl that this 
more realistic pair binding is experimentally observable 
in some cases. 



C. Cooper pair bound states from the 
Bethe-Salpeter equation 

Neglecting the phonon dynamics in electron scattering 
is actually not justified in band-insulators or low-density 
metals. A typical momentum transfer can be arbitrarily 
small in the absence of a Fermi surface, and thus smaller 
than the energy transfer from one electron to another in 
a scattering event. Given that the speed of sound is small 
on the scale of typical electron velocities, retardation ef- 
fects caused by the exchange of phonons could be quite 
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FIG. 3: The diagrammatic Dyson equation for the two- 
particle Green's function Q. 



large. We can take them into account only by the proper 
field-theoretical treatment of scattering. Therefore, we 
will solve here the Bethe-Salpeter equation for the 2DES 
to find the Cooper pair bound state. This is a relatively 
complex task, so we will make certain approximations to 
reach the important qualitative conclusions quickly. We 
will find that the retardation and quantum effects reduce 
but not eliminate the short-range attractive potential re- 
sponsible for the bound-state. By considering the worst 
case scenarios, we will argue that the effective short-range 
attraction is weakened from the one used in the previous 
section by a factor of two or so. 

The Bethe-Salpeter equation is obtained from the 
Dyson equation for the two-particle Green's function, 
whose diagrammatic representation is shown in the FigJ3] 



G = Go + G rg 



(23) 



Technical details behind its derivation and properties are 
reviewed in the appendix The full two-body Green's 
function G, the bare Green's function Go and the interac- 
tion vertex V are represented by matrices whose rows and 
columns are indexed by (P M ,p M ), where P M = (P, Q) = 
Pin + P2n is the total and p^ = (p,w) = (pi M - f>2 M )/2 
the relative momentum/energy of the two electrons in a 
Cooper pair. If there is a two-particle bound state with 
total energy £p, then the Green's function has a pole as 
a function of 57 and can be separated into a coherent and 
incoherent part: 



G(P fi ,P fl ;P' fV P , l J = 
x (2^) 3 <5(P 



F(P;p u )^(P';p' 



n-£ P + i0+ 



(24) 



p' 



The total momentum and energy are conserved, hence 
the 6(Pu — P'fj) factor. This conservation law also guar- 
anties that the function J^P,^) = F cm (V)F(p ll ) from 
the coherent term is a product of two functions associ- 
ated with the center-of-mass (.F cm ) and relative (F) mo- 
tion of the Cooper pair's electrons. We will isolate F(p IJr ) 
because it is related to the Fourier transform ip(p) of the 
stationary bound-state wavefunction ip(r); 



V>(p) = J 7^ F (P^ 



(25) 



The bare two-body Green's function is given by the prod- 



uct of two single-particle propagators ([3]) : 



GuiP ll .p ll :P' ll .p' ll ) = iG[ ^ 



Pp. G 



Pn 



Pi* 



x(2^ 3 ,5(P A1 -P;)(2 7 r) 3 ^-^) 



(26) 



and conserves the momenta and energies of both par- 
ticles. The single-particle propagators should generally 
contain all self-energy corrections, but in this case we 
will ignore them because the weakness of the electron- 
phonon coupling justifies working at the lowest order of 
perturbation theory. This also means that the vertex 
part r need not contain more than one phonon line and 
two electron-phonon vertices. If the phonons were able 
to propagate only in the two dimensions of the quantum 
well, their vertex part 



r 2Z5 (P M ,^;P;,^) = g 2 D{p^ -^)(27r) 3 <J(P M , „ 



P') 
(27) 

would involve only the (2+l)D momentum and energy 
transfer p^ — p'^ between the two scattering electrons. 
However, phonons can generally propagate in all three 
dimensions, so their propagator D(6p fl ) should take a 
(3+l)D momentum/energy transfer. The extra momen- 
tum component q z perpendicular to the quantum well is 
not conserved by the coupling of phonons to the quan- 
tum well electrons. We shall regard the quantum well as 
being very thin and limit the momentum q z by its own 
cutoff A z ~ d~ l , where d is the quantum well thickness. 
We will absorb any geometric reduction of the electron- 
phonon coupling into g<j. The in-plane momentum com- 
ponents are bounded by the cutoff A ~ a -1 . Phonons 
with sufficiently large wavelength along the z-direction 
will couple to the electrons in the quantum well in es- 
sentially a momentum-independent manner. Therefore, 
we can approximate the vertex part of three-dimensional 
phonons as: 



?MPu,Pu\ P'M = So / ^D(p p' + iq z 



x (2n)*6(P^ P' ) 



LU — Uj) 

(28) 



In order to simplify the initial discussion, we will carry 
out most derivations assuming that the phonons have a 
purely two-dimensional dispersion captured by the vertex 
part r = T2D ■ It will be straight-forward to switch later 
to the three-dimensional phonon dispersion. The Bethe- 
Salpeter equation for the relative motion P(p M ) of two 
electrons in a Cooper pair is obtained from (|23l) and (p4|) 
by taking the limit — > £ p and neglecting the terms 
that do not diverge on the approach to the bound-state 
pole: 



F = QqTF 



(29) 



In this equation we canceled out the common factors that 
depend only on P and now index all vectors and matrices 
by the relative (2+l)D momenta p M . The expanded form 



9 



of this equation involving two-dimensional phonons is: 



F( Pfi ) = ig 2 G 



P„ 



G 



P„ 



(30) 



(2tt) 3 



D(p-p',uj-uj')F(p') 



If one could drop the frequency dependence of the phonon 
propagator D, and thus neglect energy transfers among 
electrons in scattering events, then the existing integral 
over uj' of F(p',uj') on the right-hand side would pro- 
duce precisely the desired stationary wavefunction ip(p) 
according to ([25]). Integrating out the whole equation 
over uj would reduce the left-hand side to ip(p) as well, 
so one would be left with an integral form of a station- 
ary Schrodinger equation. This ui integral would pick 
a pole from one of the single-particle Green's functions 
on the right-hand-side, and allow a straight-forward con- 
version of the above integral equation to the standard 
differential Schrodinger equation. The outcome of this 
procedure would be precisely the equation we solved in 
the previous section. However, we will have to adapt this 
procedure to the case of non-negligible energy transfers, 
where strictly speaking the ordinary Schrodinger equa- 
tion is not an accurate description of dynamics. The 
Schrodinger equation can describe only the stationary 
energy levels, but not the phonon-mediated transitions 
between them which are captured by the Bethe-Salpeter 
equation. The only bound state of a Cooper pair that 
remains stable and stationary despite energy transfers is 
the ground state. Therefore, we may sensibly reduce the 
Bethe-Salpeter equation to an effective Schrodinger equa- 
tion only for the ground state. This is sufficient for the 
purpose of estimating the binding energy. The reduction 
will incorporate all corrections to the ground state due 
to the zero-point quantum fluctuations of phonons (we 



will assume that there is no external source of phonons, 
and work at zero temperature). 

We are interested in the Cooper pair ground state at 
rest. Therefore, let us substitute the propagators (@| and 
P = in (El: 



ig- 



Ml 4. 

2m ~ 

(2tt) 3 (w - uj') 2 - u 2 \p - p' 



(31) 



u 2 \p 



,'|2 



i0" 



-F{P'») 



The electron chemical potential \i from in (j4| is neg- 
ative in the band-insulating state, meaning that it takes 
a finite energy to create a quasiparticle excitation. 
We have introduced here the binding (potential) energy 
e = £p — |P| 2 / (4m) — 2|/i| < of a Cooper pair with total 
energy fl = £p. 

The structure of ([31]) suggests that we should seek a 
solution of the following form: 



F(p,u) = 



JpP 

2m 



i0+ 



IpI 2 

2m 



(32) 

where $(p, w) is a smooth function of ui without singu- 
larities in the upper complex half-plane. The numerator 
is chosen to cancel out the residue of the frequency in- 
tegral in ([25]) and establish a simple relationship to the 
wavefunction: 



iHp) = *(p,!-St 

1 2 2m 



(33) 



We will thus set uj = e/2 
above form of F(p, uj) in 



p| 2 /2m and substitute the 
D: 



IPT 

m 



V>(p) 



d p 



I' 



d> 



u 2 \p- P '\ 2 



*( £ -^f) $ (P'^') 



(-ML _i_ Jp 

^ 2m ' 2j 



(2tt) 3 (w - uj') 2 - u 2 \p - p'\ 2 + i0+ 

<Mp') 



Ip'P 



+ uj' + i0+ 



Ip'I 2 



i0+ 



« 2 lp-p'l 2 



l^ 2 



a 

2 



2„/ 



dV 



U 2 |p- P '| 2 + j0+ 



(34) 



, _ mi _ ML _ ip _ p/i + i0 + 

2m 2 m \" " I 



^ + ^+u|p-p' 

2m 2m 1^ 1 



i0+ 



* P' 



2 2m 



wp-P 



We have carried out the integral over uj' on the right- 
hand-side of ([3"Tjl and obtained two terms. One term is 
generated by the particle pole and contains the wave- 
function ip(p'), while the other term is generated by 



the phonon pole and still contains the unknown func- 
tion $(p',w') at the frequency uj' that depends on the 
unknown energy e. Dealing with the latter is quite dif- 
ficult, so we will introduce our first approximation. The 
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"phonon-pole" integral in the last line of (|34|) is domi- 
nated by the values of p' at which: 



2m, 



IPT 

2m 



up 







(35) 



This is the only way for the integrand to diverge, given 
that we expect no singularities in $(p', u/) and e < 0. We 
will substitute the above special value of |p'| everywhere 
in the "phonon-pole" integral, except for its diverging 
denominator. This yields a very simplified approximate 



contribution of the phonon pole: 



2„/ 



d Z p 



up 



>(p') (36) 



2 j -ff+e+«ip-p' 

whose main benefit is that the unknown $(p',u/) is 
traded for the wavefunction tp(p') that we wish to de- 
termine. This approximation reduces the Bethe-Salpeter 
equation to: 



IPT 

m 



^(p) « g 2 



d p 



u 2 |p- 



VI 2 



V 2m 2m 



+ 



up - 



« 2 |p - 



VI2 



2 m 



IP'I 2 

2m 



u|p-p' 



(37) 



At this point we have an integral form of a Schrodinger 
equation for the Cooper pair bound state, written in 
the momentum representation. We can easily convert it 
to the position representation by performing the inverse 
Fourier transform (note the appearance of the reduced 
mass to/2 in the kinetic energy): 

V>(r) + g 2 I dV V{r, r')V'(r') = eip{r) . (38) 



The only price to pay is that the potential V(r, r') is non- 
local due to the retardation effects. We can view V(r, r') 
as a non-diagonal operator in the position representation. 
Its components are: 



V(r,r') 



^^e^ P V^ (p)p0) (39) 



where V(p,p') is the expression from the square brack- 
ets of (|3T|) . Upon a closer inspection of V"(p, p') it can 
be easily seen that the potential energy operator V is 
not Hermitian. The eigenvalues e of the ensuing non- 
Hermitian Hamiltonian will generally end up being com- 
plex. This indicates that Cooper pair states have a finite 
lifetime of the order of 1/Im(e). Indeed, Cooper pairs can 
emit phonons and relax to lower energy states. Only the 
ground-state is stable and its eigenvalue e is real. Being 
interested only in the ground-state, we will symmetrize 
the potential energy operator 



l r 



y s (r,r') = - U(r,r') + ^*(r',r) 



(40) 



and hence eliminate its non-Hermitian parts. This will 
have no effect on the Hamiltonian projection to the 
ground-state, and the ground-state energy e. The real- 



space symmetrization is equivalent to: 
1 U 



V a {p,p') 



2 

; 2 |<5p| 2 



V(p,p' 



V*(p',p) 
u\Sp 



(41) 



SE 2 



u 



■\5p\ 



u\6p\ 



l \Sp[ 



2 SE 2 



*\6p\ 



5E + u\5p\ -5E + u\5p\ 
P SE 2 



2 5E 2 



'■\Sp\ 



where P stands for the principal part, and we defined 
5p = p p' and 5E = (|p| 2 — |p'| 2 )/2m for brevity. 
The symmetrized potential has two components. The 
first component —1/2 is momentum-independent, so it 
corresponds to a purely local attractive contact potential 
in real space. We have shown in the previous section 
that such a potential produces a bound-state. The second 
component is momentum-dependent and non-local in real 
space. It exists because we included finite momentum 
transfers (SE ^ 0). Since it has a repulsive character for 
some momenta, we must check if it can weaken or destroy 
the bound-state. 

After symmetrization we are guarantied to get real 
eigenvalues e. However, the non-local Schrodinger equa- 
tion is still too difficult to solve. We will now make our 
second "self-consistent" approximation, which will allow 
us to verify the survival of the bound-state in the pres- 
ence of the non-local component of the potential. We 
will start with a hypothesis that the ground-state of two 
electrons is indeed a bound state, and use it to construct 
an effective purely local potential that approximates the 
realistic one. This approximation will overestimate the 
non-local component of the full potential, but accurately 
include its local component that by itself would support 
a bound state. Therefore, if this effective local potential 
has a bound state, so does the full one and the starting 
hypothesis is justified. 

The assumed bound state has zero angular momen- 
tum, so its wavefunction amplitude is largest at the zero 
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separation |r| =0 between two electrons and gradually 
decreases as |r| grows. Consequently, replacing i/j(r') in 
(|38f by ij)(r) would overestimate the effect of the non- 
local component of the full potential at |r| < a. If the 
non-local part is repulsive, this overestimate would yield 
an upper bound for the binding energy e < 0. On the 
other hand, if the non-local part is attractive at short 
distances, the overestimate would produce a deeper bind- 
ing energy than the real one, not particularly useful for 
quantitative purposes. In either case, however, a bound- 
state is guarantied to exist if the overestimate supports 
it. Note that this procedure would clearly not modify the 
purely local part of the potential, which is short-ranged 
and attractive according to (|4Tj) . 

We will actually self-consistently find a bound state 
only when electron interactions are mediated by phonons 
that disperse in all three dimensions. Purely two- 
dimensional phonons that we currently pursue produce 
a repulsive effective local potential at short distances, 
which wc will derive first. Replacing ip(r') by iP(t) in 
the symmetrized yields the following effective local 
potential 



U(r)= / dW s (r,r') = 



d 2 p 



(27T) 



2 e*P r ^(p,0) (42) 



d p 



JeI! 

2m 



IpI 2 

2 m 



u 2 |p| 2 



and a standard Schrodinger equation 
V 2 

ip(r) + g 2 U(rU(r) = eipM . (43) 

m 

The local potential (|42p created by the two-dimensional 
phonons can be easily calculated: 



u 2 \p\ 2 



l €) -" 2 ipi 2 



l [dp P 

n \1raul 



(mu)' 



Y {2mu\v\) , 



(44) 



where J n (z) and Y n (z) are the Bessel functions of the 
first and second kind respectively. This potential is re- 
pulsive and logarithmically divergent at |r| —> 0. It is 
also oscillatory at larger distances, alternating between 
positive and negative regions. A potential of this kind 
could support bound-states only if the typical electron 
momentum scales were of the order of mu. This is, how- 
ever, not the case. The sound velocity u is much smaller 
than the typical electron velocity. 

The coupling of the quantum well electrons to the bulk 
phonons of the superconducting material is captured by 



the vertex (j28|) . and this gives rise to a slightly different 
Bcthe-Salpeter equation ([50]) : 



P„ 



(45) 



(2tt) 3 2tt 



D(p-p' + zq z ,u J -oj')F(p'J 



The electron-phonon coupling g Q has different engineer- 
ing dimensions than g to match the changed dimensional- 
ity of the integral. We can immediately sec that the pro- 
cedure of deriving an effective Schrodinger equation can 
be translated directly to the case of three-dimensional 
phonons merely by changing the phonon dispersion from 
«|p — p'| to u\p — p' + zq z \ and integrating out q z in 
addition to p'. Adapting (l42l in this manner yields: 



U 3D (r) 



d 2 p dq z e lpr 



(2tt) 2 2tt 2 



(mix 

\2m J 



IpI 2 

2ni 



u 2 \p + zq z \ 2 



2na 2 d 



(46) 



where 



% D (v) = 



dp dq z e lpr 



fJpl!V 

I 2m J 



(2tt) 2 2t: 2 f \ p \ 



\ 2m J 



ii 2 |p| 2 — u 2 q 2 



2mu 

Au 2 J 2tt \2m 
o 

3(mu) 3 



Jo(p\r\) 



1 - (-£-) 

V 2mu I 



1( . iF 2 I |;l,3;-(m«|r|) a 



(47) 



Therefore, the bulk phonons create the core attractive 
potential at |r| < a ~ A -1 that we analyzed before, and 
an additional oscillatory potential U' 3D {r) given by the 
generalized hypergeometric function 1F2 whose asymp- 
totic behavior is: 



1F2 [ 3! M; 



1 



-x 2 

6 



35 4 
192 A 



U.x-', 



sin(2x) 




The extra potential U^ D is actually attractive at short 
distances, but recall that we have overestimated its con- 
tribution. At any rate, the overall scale (mu) 3 of this 
potential is small in comparison to that of the attrac- 
tive core, which is determined by the momentum cut-offs. 
The most significant correction from the naive analysis in 
the previous section is that the amplitude of the attrac- 
tive core potential is reduced by a half due to the phonon 
dynamics (recall that A ~ A /d = g 2 /d). The Cooper 
pair bound-states survive and have the same binding en- 
ergy by the order of magnitude that we calculated before. 

So far we have treated the 2DES as being electroni- 
cally isolated from its environment at all energy scales. 
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However, in reality we hope to engineer the electron- 
phonon coupling in the 2DES by the proximity effect. 
This means that the 2DES is coupled to a bulk super- 
conductor whose spectrum contains extended states at 
energies \E — fi\ > A , where Ao is the superconductor's 
pairing gap. Only those electrons in the 2DES whose 
in-plane momentum is small enough to give them energy 
\E — fi\ < Aq can experience the true two-dimensional 
dynamics. Hence, we can patch our previous results by 
introducing a different momentum cutoff, p < Aa <C A, 
which is related to A . Only the core attractive potential 
in pp]) depends on the cutoff, while such dependence of 
(|44| and (|47|) can be neglected. The patched core po- 
tential is weaker and has a larger range l Ao ~ a 
than the one written in (|4"r?)) . This effect is due to the 
ability of electrons to tunnel across the superconductor- 
2DES interface and completely unrelated to the electron- 
phonon coupling A. Thus, it will plague even strongly 
coupled superconductor interfaces. We will consider the 
naive formal A — > oo limit to isolate this issue. The en- 
suing Cooper pair binding energy scale Eq at A — > oo 
must be much smaller than the naively obtained scale 
1/roo 2 in the section Hi Bl Note that the binding energy 
reduction to Eg has little to do with the increased range 
of the core potential in the weak-coupling limit. The 
reduction is mainly the consequence of the weaker at- 
traction strength. Thus, we can see without any further 
calculation that Eq cannot be (much) larger than Ao be- 
cause this would not allow both the fermionic quasiparti- 
cle and Cooper pair two-dimensional bands to fit within 
the bulk superconductor's pairing gap. Some (small) en- 
hancement of Eq can perhaps occur due to the in-plane 
momentum conservation. The energy scale of retardation 
corrections is still much smaller than Eq ~ Ao. 

With all this in mind, we conclude that the bulk lon- 
gitudinal phonons of the superconducting material are 
significant mediators of the interactions among electrons 
in the quantum well. It is possible, though not neces- 
sary, that a branch of two-dimensional interface phonons 
exists in the heterostructure. Even though they create 
a logarithmically divergent repulsive potential between 
electrons at short distances, this potential is unlikely to 
jeopardize the existance of bound-state Cooper pairs due 
to its overall small scale and slow divergence. Note that 
this divergence is cut off below |r| < a and thus ul- 
timately innocuous, just like the much more dramatic 
divergence of the Coulomb repulsion which we consider 
next. 



D. Coulomb interactions 

Our analysis has so far ignored the Coulomb interac- 
tions between electrons. Since they are the main adver- 
sary of the phonon-mediated pairing, we will investigate 
here whether they destroy the Cooper pair bound-state in 
two dimensions. The answer to this question is ultimately 
microscopic and system-dependent. However, we will ar- 



gue that the condition for the bound-state survival in the 
2DES is universally close to the condition for pairing in 
the superconducting material that creates the proximity 
effect, assuming that the 2DES thickness d is sufficiently 
small. In other words, an insulating 2DES close to its su- 
perconducting transition will feature low energy bound- 
state Cooper pair excitations if the coupling between its 
electrons and the proximate superconductor's phonons is 
not weakened too much by the interface geometry. 

Let us first analyze the condition for pairing in con- 
ventional BCS superconductors. The phonon-mediated 
interaction between electrons can be accurately captured 
by the Bethe-Salpeter equation with the phonon prop- 
agator in ([3]). Wc may neglect phonon-mediated en- 
ergy transfers between colliding electrons in weakly cou- 
pled superconductors (wd <C Ef), and approximate the 
phonon propagator D(p^) with the Fourier transform 
^bcs(p) °f t ne effective electron-electron interaction po- 
tential L/g CS (r) written in ([5]): 

^bcs(p) = . (49) 

This expression is valid only below the cutoff momentum 
scale A ~ a -1 , where a is the inter-atomic spacing. We 
can view Ao as being roughly constant at p < A and 
zero at p > A. Interactions at larger momentum scales 
are negligible because they can be generated only by the 
electron coupling to higher phonon bands. Interacting 
electrons near the Fermi surface exchange momenta of 
the order of pf and smaller, which lie sufficiently below A 
to allow treating Ao as a constant. The screened Coulomb 
interaction strength in momentum space is: 

^cs(P) = ^. (50) 

where l s = p^ 1 is the screening length and we use the 
Gaussian units. Through l s we include the effect of 
all electrons on the net repulsive interaction between 
any two electrons. This interaction is not subject to 
any momentum cutoff since the corresponding real-space 
Coulomb potential E^bcsM = e V r snou ld hold at arbi- 
trarily small r C I s . We could in principle improve on 
this expression by incorporating the relativistic correc- 
tions and some elements of photon dynamics, but this 
is unnecessary because the typical electron velocities are 
sufficiently smaller than the speed of light. 

The universal prerequisite for pairing is that the to- 
tal electron-electron interaction be attractive at impor- 
tant momentum transfers which occur in most scattering 
events. This condition is 

^BCs(p) + t / Bcs(p)<0 , p<Pt. (51) 

in Fermi liquids. Let us take the superconducting Nio- 
bium as an example in order to estimate the relevant 
scales. In Niobium, E f = 5.32 eV, p f = 1.24 x 1CT 24 kg • 
m/s (the effective mass m* is approximately the same as 
the free electron mass m), a = 3.3A and A = 2ixja — 
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2.01 x 10 24 kg • m/s. We can estimate the screening 



length using the Thomas-Fermi formula^: 



Ps 



(52) 



which yields p s = 1.77 x 10~ 24 kg • m/s. This is a typi- 
cal low-temperature situation in metals and conventional 
superconductors based on them, pf sa p s ss A. Knowing 
this, we can write a simple approximate form of the con- 
dition (EB: 



An > 



4ire 2 



(53) 



The electron-phonon coupling must be sufficiently strong 
to give rise to superconductivity. Since the dynamics 
in BCS superconductors involves momentum exchanges 
near or below pf, where (|50p exhibits very little mo- 
mentum dependence, we may use the effective coupling 
Aq = An — 47re 2 /p 2 in BCS formulas such as ([7]) to obtain 
the pairing gap, critical temperature, etc. This would 
include the net effect of phonon-mediated and Coulomb 
interactions on superconductivity. 

Now let us obtain the analogous condition for the exis- 
tance of bound-state Cooper pairs in an insulating quan- 
tum well that has an interface with a conventional super- 
conductor. The phonon-mediated electron-electron inter- 
action will be again represented by a static momentum- 
independent potential: 



U P (P) 



-X 



(54) 



that holds at p < A. The effective coupling A between the 
2DES electrons and the superconductor's phonons has 
different physical dimensions than the superconductor's 
An because U p (p) is a two-dimensional Fourier transform 
of a short-range potential. A crude estimate A < And -1 
involving the 2DES thickness d is based on the physical 
dimensions and the fact that the heterostructure geome- 
try plays the main role in reducing the value of A below 
the one set by the bulk electron-phonon coupling in the 
superconductor. We are neglecting the retardation ef- 
fects due to phonon dynamics, because we found in the 
previous section that they are small. 

The effective Coulomb potential in momentum space 
has different behaviors at small and large momentum 
transfers. Small momentum transfers correspond to large 
distances between electrons. Since we are considering an 
insulating quantum well, its electrons cannot give rise to 
screening. However, the superconductor's electrons can 
partially screen the quantum well electrons by creating 
"image charges" at the interface. Each electron excita- 
tion in the quantum well is effectively turned into a dipole 
whose moment points in the direction perpendicular to 
the well and has magnitude ~ ed, or ~ e(d + /b) if the 
tunneling barrier is present (we will assume no tunneling 
barrier in the following discussion). The ensuing dipolc- 
dipole interaction is repulsive: 



Ce(r) 



2 j2 



e 2 d 



r > p s 



(55) 



Note that the superconductor cannot efficiently provide 
"image charges" at electron separations r smaller than its 
screening length l s = p^ 1 . The two-dimensional Fourier 
transform of this potential is: 



^ / ^ e " Pr - 27r(ed) 2 / dr^M 



U c (p) ~ eV / d^r ; 



27r(ed) : 



1*2 | -^}^ 



1 1 




Ps -P 



^2n(ed) z ^ Ps -p+^-j , p<p s . (56) 

The last line is obtained from the Taylor expansion of 
the generalized hypergeometric function 1F2 in the limit 
where this expression applies. At momentum transfers 
sufficiently above p s we must go back to the unscreened 
Coulomb potential 



U c (r) ~ - 
r 

whose Fourier transform is 



r < p s 



tt l ^ 2ire2 p 1 1 3 f P " 



(57) 



(58) 



2^e 2 
P 



(small oscillations) , p > p s 



Using the previously developed picture of the Cooper 
pair bound-state in the absence of Coulomb interactions, 
we would conclude that the main momentum transfers 
between the two electrons of a bound-state Cooper pair 
occur just below p ~ a" 1 ~ A in a relatively broad range 
of momenta (down to k ~ £ -1 )- The effective electron 
repulsion is approximately given by (|56[) in this range of 
momentum transfers, and this conclusion is not altered 
by the reduced cutoff Aa due to the interface between 
the 2DES and the bulk superconductor. Even in the 
worst case, the Coulomb interaction effects can be weak 
enough for the bound-state to survive (U p + U c < 0): 



A > 27re 2 d 2 



(59) 



We can compare this with ([53]) by converting the 2D cou- 
pling A to the equivalent 3D coupling Ac? which is quanti- 
tatively comparable with the superconductor's electron- 
phonon coupling An or slightly smaller when d is suffi- 
ciently small. We will also use the fact that pf sa p s sa A, 
and assume that the quantum well thickness d is larger 
but not much larger than the typical inter-atomic separa- 
tion ~ pj 1 of the superconductor. Then, we can rewrite 
([59"]1 as: 



27re 2 , 
\d > — —(pfdy 

Pf 



(60) 
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Since Ac? is of the order of Ao , we see that the conditions 
(f53"| and ([50]) are not too far from each other. If the geom- 
etry of the coupling between the electrons in the quantum 
well and the superconductor's phonons does not make A 
too small, the existance of pairing in the superconductor 
can imply the existance of bound-state Cooper pairs in 
the quantum well despite the Coulomb repulsion. 

It should be noted that the net interaction potential re- 
mains dominated by the Coulomb repulsion at extremely 
short distances between electrons. This might appear to 
pose a problem given that a bound state in two dimen- 
sions is guaranteed only when the potential is attrac- 
tive at short distances. However, such short distances 
in our system can be probed only by scattering events 
whose momentum transfers lie above the cutoff. This 
in fact corresponds to electron scattering between differ- 
ent bands, which is suppressed at low temperatures. The 
Cooper pair bound-state involves low-energy band quasi- 
particles, and is shaped at energies and momenta that lie 
below the cutoff. 

We can now apply the above analysis to make some ex- 
tremely crude estimates of the critical temperature T 2D 
for pairing and superconductivity in the 2DES, knowing 
the critical temperature T c of the proximate supercon- 
ductor material. Assuming weak pairing in the super- 
conductor, we can extract its effective electron-phonon 
coupling from the formula ([7} using^ Aq ~ 1.76T C : 



By combining the last four formulas we obtain: 



a;, 



2tt 2 



m Pf lo s(o%) 



(61) 



We will treat the value obtained from this formula as 
empirical, which means that it incorporates both the true 
electron-phonon coupling Ao and the Coulomb repulsions 
from dST 



A ~ Aq — 



47re 2 



(62) 



The intrinsic proximity-generated coupling between the 
2DES electrons and superconductor's phonons is of the 
order of A ~ Ao/rf, but the effective coupling A' that 
determines the critical temperature T 2D is also reduced 
by the Coulomb repulsions according to (|6"U|) : 



A' = ^ 



27re 2 



{Pidf 



(63) 



Finally the critical temperature in the 2DES is related to 
the Cooper pair binding energy (l22l) : 



T: 



2D 



Eq exp 



47T 

mA' 



(64) 



r c 2D ~^ exp 



Eq exp 



2p f d 



*/Mis&) + a ^[i-*(w9 I 



(65) 



Clearly, T 2D is maximized in the d — > limit, but will 
not exceed T c because Eq ~ Ao. There is also a critical 
2DES thickness d at which T 2D -> 0: 



1 

Pi 



me 2 log ^ 



76T< 



(66) 



If d > G? max , then the Coulomb interactions win in the 
2DES and Cooper pairs cannot be formed. 

The Niobium parameters pf — 1.24 x 10~ 24 kg • m/s, 
T c w 10 K « 0.86 meV, uj d « 275/2tt K w 3.77 meV 
yield <i max ~ 4.89 A. This may be too small, but is not 
too far from the more realistic values d ~ 10 — 100 A. 
One must keep in mind that this estimate is extremely 
crude and does not by itself rule out Niobium as a ma- 
terial that can produce an observable proximity effect. 
An ab-initio microscopic calculation would be able to 
produce a more accurate estimate, but only experiments 
can really verify practicality. The present crude analy- 
sis is more useful for comparison purposes. Consider, for 
example, the best conventional superconductor MgB2- 
Its relevant parameters^ - — are T c sa 40 K sa 3.44 meV, 
E f ~ 0.25 eV, o; D « 1100/2tt K - 15.09 meV. This 
material with multiple and anisotropic bands is not ac- 
curately described by the simple weak-coupling BCS the- 
ory, but we will nevertheless use the BCS relationships 
to make ballpark estimates. Taking the worst-case sce- 
nario m* « m gives the upper bound on the Fermi mo- 
mentum pf ~ y / 2mE{ w 2.7 x 10 -25 kg • m/s and the 
worst conditions for superconductivity in the 2DES. Now, 
d max ~ 29.49 A. A quantum well of thickness d = 20 A 
would yield a critical temperature T 2D ~ 0.83-Eo, which 
is close to the superconductor's T c ~ 40 K given that 
E Q - T c . In fact, the dependence T 2D (d) is very weak 
until d ~ d max because the factor l/mpfe 2 d 2 lax ~ 0.09 
that appears in (|65l) is very small for MgB2. Note that 
even if this estimate were off by a factor of 100, the en- 
suing pairing temperature scale ~ 400 mK would be de- 
tectable. Therefore, while all these approximations can- 
not be taken seriously for quantitative purposes, they in- 
dicate by orders of magnitude that the proximity-induced 
superconductivity via the phonon mechanism can be 
achieved as a matter of principle in realistic quantum 
well heterostructures. 
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III. COOPER PAIR TUNNELING 

Here we derive the effective action of a two-dimensional 
electron system (2DES) placed in contact with a bulk 
superconductor (SC) by focusing entirely on the virtual 
Cooper pair tunneling from the SC to the 2DES. We will 
not pay attention to phonons in this section, but their 
contribution to the effective action can be easily included 
via an extra two-body potential that we already discussed 
in detail. Our present goals are to find how the Cooper 
pair dynamics renormalizes the electron dynamics in the 
TI, and demonstrate that even the tunneling mechanism 
produces short-range attractive interactions which have 
a spin-independent component capable of forming triplet 
pairs in the 2DES. 



A. Singlet proximity effect 

Many aspects of the interaction proximity effect can be 
derived from the continuum limit path integral in which 
we model the SC as a material with instantaneous short- 
range attractive interactions among electrons, and the 
2DES as a multi-orbital band-insulator without any in- 
trinsic interactions. We will use t to denote imaginary 
time in this entire section, and r to label the 2DES 
electron's orbital degrees of freedom. The SC-2DES in- 
terface in the xy plane is modeled simply by orbital- 
dependent electron hopping t T between the SC and 2DES 
surfaces. The simplest imaginary-time action of this kind 
is S = S p + S g : 



dt d 3 r 



z<0 



fl 



D l - 

2m 



Use f a 



dt d 2 r 



2DES 



^ \D + H U 



J dzY,U{z)[^ Ta l^ +z J a (z) + h.c) 



z<0 



S e = jdtd 3 r-F^F^ , 



(67) 



where the SC's electron fields f a are defined in the semi- 
infinite space z < 0, while the TI's electron fields ip Ta 
are two-dimensional with spin a and orbital r degrees of 
freedom. The gauge field dynamics is governed by 
S g via the electromagnetic field tensor in the entire 
space, while the coupling of matter to the gauge field is 
embedded in: 



Dp = dp - ieA^ 



I12 = exp 



/ 



V 



dx^A^ 



(68) 

Specifying the bare Hamiltonian H of the 2DES's elec- 
trons is not essential for the present purposes. The SC 



and 2DES are coupled via the hopping term t T (z), whose 
magnitude depends on the hopping distance. The pos- 
sible spatial separation of the 2DES's orbitals along the 
z-axis translates into the r and z dependence of hop- 
ping, while we neglect hopping processes in which the 
electron's x and y coordinates change (without missing 
any essential physics). 

We wish to integrate out the SC's fields and obtain 
the effective action for the 2DES alone. To that end, we 
first apply the Hubbard- Stratonovich transformation and 
decouple the interaction in the SC using a singlet Cooper 
pair field <f>. Then, we use the Bogoliubov transformation 
to group together all terms quadratic in electron fields: 



z<0 

where 

"HBdG 



■S',, = J dt d 3 r i^ 2 + ( ^ f i ) H B dG f f f \ J + J 
<o L V 4- / J TI 



Y Do + ^ + /i sc 

The Gaussian integration of f a leads to: 
" 1 



(69) 
(70) 



Sf = / dt d 3 r 



z<0 



jj\<t>\ 2 ~ trlog"H B dG 



(71) 



d 2 rdt 



2DES 



^ (D + Ho) V 



z<0 




We have redefined tl> TrT (r, z) — > tp Ta (r)l r r+z± in all three- 
dimensional integrals in order to maintain gauge invari- 
ance. The <f> field implicitly contained in "HedG is three- 
dimensional, but mediates interactions between the two- 
dimensional 2DES's electrons. We can invert the matrix 
%BdG by representing it as: 



Hbag = Do + T + $ 



(72) 



Di 



where a a are Pauli matrices in the Nambu spinor space. 
The matrix $ now contains all spatial variations brought 
about by the fluctuating field 0, and we can express 7ig] G 
as a gradient expansion: 

^BdG = (Do + T + S)" 1 (D — T— $r x (D - T - $) 
= (Dl - T 2 - $ 2 + Stpy 1 (D - T - *) , (73) 

= f E P") (D 2 T 2 $ 2 ) 1 (Do T *) , 

\n=0 / 
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where 



electrons: 



V = - (D 2 - T 2 - $ 2 ) 1 <5$ 



<5$ = (D - T - $) (D + T + $) 
= [A),T + $]-{T,$} . 



(74) 



(D 2 -T 2 -<f> 2 ) 



li 4> = const then = 0, so the operator T> gathers only 
spatial and imaginary time (gauge-covariant) derivatives 
acting on (j>. Hence the effect of V can be neglected deep 
in the superconducting phase. Keeping only the n = 
term amounts to the mean-field approximation. We can 
also consider the effect of D 2 , and kinetic energy T 2 small 
next to $ 2 = 1 4>\ 2 (dominated by the SC's pairing gap). 
If we substitute H^dG mto we reveal the generated 
pairing interactions in the 2DES at the lowest order of 
the gradient expansion: 



Sf = I dt d 3 r 



z<0 



tr log H B dG 



(75) 



+ J Atd 2 rl^(D + H )i, 

2DES ^ 



h.c.) H 



z<0 



The dots denote various terms involving the derivatives 
of cj>, and renormalizations of the bare 2DES Hamiltonian 
H . 

In the last step we integrate out <j> and thus reverse 
the original Hubbard-Stratonovich transformation. The 
dynamics of <j> is shaped by: 



Sa 



did'V 



\cj>\ 2 - tr log "HBdG 



(76) 



z<0 



dt d A r 



z<0 



0t (do - 2ieA )(j> + —\{di- 2ieA t 



t\4>\ 2 + uW 



The general form of this action and the charge of <fi are 
restricted by symmetries, but calculating the concrete 
values of the couplings such as t and u is not worth the 
needed effort. We will be satisfied by the knowledge that 
t > carves a potential minimum at a finite \<j>\ = \<j>o\ 
due to SC's superconductivity. The value of \4>o\ (in the 
units of energy used here) is given by the pairing gap in 
the SC material. Assuming that the SC material is not 
close to its phase transition to the normal state, \<po\ is 
reasonably large and only small fluctuations 54> = 4> — 
<po have a significant dynamical effect. Integrating out 
Sep produces the following effective action of the 2DES's 



ceff _ 



dtd 2 r 



2DES 



^ (Do + H') i> 



(77) 



z<0 




X ^ ^Tl ^T2 ^T3 ^T4 

I0o| 4 



The operator H' contains the renormalized orbital- 
mixing kinetic energy and the dots represent the ef- 
fect of the neglected higher derivative terms and large- 
amplitude 5(f) fluctuations. The correction to Ho in H' 
is small compared to the intrinsic 2DES Hamiltonian Ho 
when the SC's gap \4>o \ ~ A is large. 

The generated interaction in the 2DES depends on the 
renormalized interaction coupling U', not the bare one 
U. The renormalization comes from (1761) and pertains to 
the stiffness of order parameter fluctuations in the vicin- 
ity of <f>o. We see that the generated interaction is at- 
tractive, but its strength decreases when the SC's gap 
|^o | grows. This may look counter-intuitive for a prox- 
imity effect, but is in fact natural having in mind the 
dynamical origin of the considered interaction. After all 
approximations, we have in fact obtained a perturbative 
result in which the SC-2DES hopping is a small energy 
scale in comparison to the SC's gap \<po\- The interaction 
is induced by virtual tunneling of a Cooper pair from the 
SC to the 2DES and back, which involves four electron 
hoppings and corresponds to the fourth power of hopping 
constants t T . The tunneling occurs within a thin layer of 
the SC near the interface, where t T (z) is not too small. 
Temporarily removing an electron from the SC costs en- 
ergy \4>o\, so if the SC's superconductivity is robust its 
Cooper pairs find it hard to tunnel and mediate inter- 
actions in the 2DES. Recall, however, that interactions 
between the 2DES's electrons mediated by the nearby 
SC's phonons are still nearly as strong as inside the SC, 
and thus dominate when <f)o is large. 

This derivation shows that electron tunneling from the 
SC generates attractive interactions in the 2DES that can 
give rise to both inter- and intra-orbital pairing in the ab- 
sence of a symmetry to protect against orbital mixing. In 
addition to interactions, we obtained a direct order pa- 
rameter coupling proportional to ^Qip'tp^ + h.c, which 
is usually considered in the superconducting proximity 
effect literature. The latter explicitly breaks the global 
U(l) symmetry and produces pairing correlations in all 
circumstances, so we can say that the 2DES always hosts 
a singlet condensate as a result of the proximity effect. 
However, this effect is not equally significant in all cir- 
cumstances. The actual magnitude of the order parame- 
ter \ipo \ in the 2DES also depends on the spectrum of the 
renormalized Hamiltonian H' . If this spectrum contains 
gapless modes, as in the case of a helical Fermi liquid 
formed on the surface of a strong topological insulator, 



17 



then \ipo\ can be sizable as a result of the direct order pa- 
rameter coupling. The attractive interaction would also 
immediately create a pairing instability in these circum- 
stances. The net \ipo\ can be calculated, for example, by 
applying the mean-field approximation to (1771) . On the 
other hand, if the electronic 2DES spectrum is gapped, 
such as in a topological insulator quantum well, then |<po| 
can be very small. One way to see this is to construct 
from (|77[) an effective action of fully gapped Cooper pairs 
tp ~ ipfip± m the 2DES (which could exist as excitations 
below the single-electron gap) : 



eff _ / j7j2„ 
2DES 



1 2 

ip' (d - 2ieA )ip + - 2ieAi)tp\ 



-aifoip 1 + h.c.) + t'\<p\ 2 + u'\(p\ 4 + 



(78) 



The gap implies t' > 0, so the saddle-point approxi- 
mation estimates \ipo\ ~ a\(j>o\/2t'. Recall that the pa- 
rameter a is proportional to the square of the SC-2DES 
tunneling constant, and inversely proportional to \4>o\ 2 . 
Therefore, robust superconductivity in the SC material, 
weak tunneling and a large quantum well bandgap all 
conspire to make \ipo\ small. It cannot ever completely 
vanish, however. In order to make the quantum well's 
superconductivity \(po\ large, it is necessary to turn t' 
negative. This can indeed be accomplished in a practical 
heterostructure device by exposing the quantum well to 
a biased gate across an insulating medium. An applied 
gate voltage V g can pull electrons from the SC into the 
2DES, which is captured in the effective action by the 
increase of the 2DES's chemical potential fj, oc V g hid- 
den in Hq and H'. Since Cooper pairs are charged, this 
chemical potential directly appears in t' = t' Q — 2/i. This 
is how one can tune a large \ipo\ ~ y/\t'\/2u'. We have 
emphasized only this effect in our recent analysis of the 
topological insulator quantum wells placed in proximity 
to a superconducto r 38 ! 39 . 

In the last stage of the present derivation, we wish to 
affirm that the generated attractive interactions in the 
2DES are short-ranged. For that purpose we must draw 
a distinction between the proximity effect created by a 
superconductor and a (hypothetical) superfluid. If the 
SC material were a superfluid, its matter field would not 
be coupled to the electromagnetic gauge field . Then, 
the SC's Cooper pair action (1751) with A^ set to zero 
would contain the dynamics of gapless Goldstone modes, 
and integrating them out in the following step would 
generate an algebraic long-range (Coulomb-like) interac- 
tion among the 2DES's electrons. Goldstone modes are 
gapped out by the Anderson-Higgs mechanism in real su- 
perconductors. This can be easily seen by fixing a gauge 
to pass all local phase fluctuations of the matter fields 
such as 4> on to the "longitudinal" modes of the gauge 
field Ap, via a gauge transformation: 



The action of the fluctuations 8<f> = <fi — <po derived from 
(fTo]) acquires a term proportional to e 2 \(j)o\ 2 (A ll A^) in 
this gauge when the ground-state of the SC material is 
superconducting {\4>o\ ^ 0), so all excitations in the SC, 
now described by A^, are clearly gapped. Integrating out 
Scj), therefore, produces short-range (Yukawa-like) attrac- 
tive interactions in the 2DES. 

The 2DES effective action (1771) is still written as a 



A„ 



1 



„2iA 



(79) 



gauge theory. The direct pairing term with </>oV'jV'j + ' l - c - 
looks very awkward when one recalls that the factors 
X12 from (|68[) are embedded in the electron operators 
i/>Ta(r, z). We can get rid of those X12 factors only by 
integrating out A^ and absorbing any background gauge 
field into 0o by a gauge transformation. It is convenient 
to integrate out A^ in the same gauge that we mentioned 
above in the context of the Anderson-Higgs mechanism. 
The background gauge field is then simply zero, and the 
direct order parameter coupling yields to the same form 
that is obtained in a superfluid proximity effect. Note 
that the global phase of </> cannot be eliminated by a 
gauge transformation, but its absolute physical meaning 
is dubious because it can be identified with an arbitrary 
choice of a reference phase at a single point of space-time 
(all other local phases are gauge-dependent). The physi- 
cal meaning of the order parameter coupling is that the 
local phases in the 2DES are correlated with those of the 
SC material in any gauge. 

Integrating out A^ produces additional effect due to 
the fact that electromagnetic fields are not confined to 
the SC and 2DES volumes. Photons are gapless outside 
the system we consider and mediate the usual long-range 
Coulomb interactions between electrons. The Coulomb 
forces, however, are laterally screened out in the 2DES 
since its ground state exhibits phase correlations, and 
efficiently so when \tpo\ is large as discussed earlier. Pho- 
ton screening in the semi-infinite space occupied by the 
SC material also generates an image-charge effect for 
the 2DES's electrons. Every electron in the 2DES ef- 
fectively becomes an electric dipole and interacts with 
other 2DES electrons accordingly at length-scales larger 
than the photon screening length. 



B. Triplet proximity effect 

We have seen that the proximity to a superconduc- 
tor (SC) creates via Cooper pair tunneling not only a 
direct order parameter coupling in the two-dimensional 
electron gas (2DES), but also an effective attractive in- 
teraction that qualitatively mirrors the interactions in 
the SC responsible for superconductivity. So far we have 
considered a particular spin-dependent part of this inter- 
action that leads to the pairing in singlet channels. How- 
ever, its phonon-mediated Coulomb origin introduces a 
spin-independent component to this interaction, which 
could create triplet Cooper pairs in principle if the Pauli 
exclusion could be overcome. Indeed, topological insu- 
lator quantum wells are 2DESs with orbital degrees of 
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freedom that allow spin-triplet Cooper pairs to form as 
inter-orbital singlets, and even give them energetic ad- 
vantage via the Rashba spin-orbit coupling. Here we will 
explore the ability of Cooper pair tunneling to gener- 
ate the spin-independent interaction component in the 
2DES. 

The phonon-mediated attractive interactions in the SC 
do have a small finite range established by charge screen- 
ing, which we will optimistically compare with the thick- 
ness d of the quantum well 2DES. Taking this into ac- 
count, we will find that the effective attraction inside 
the quantum well could also bind spinful inter-orbital 
triplet pairs, and we will estimate its strength. But, we 
need a more complicated model in order to capture the 
inter-orbital pairing in the quantum well. Consider the 
following action for the SC material: 



5 SC = Jdtd'r 

z<0 



(80) 



-~ J d 3 r'T/(r-r')4(r)/ CT (r)/],(r')/ CT ,(r') 

z<0 

in which we neglect the spin dependence of the attractive 
interaction V > 0, assuming its origin in the Coulomb- 
based phonon mechanism. Decoupling it by the Cooper 
channel Hubbard-Stratonovich transformation requires 
complex fields <f> aa >(r,r') that have a doubled position 
and spin dependence. In order to avoid ambiguity, one 
could temporarily reorganize the above conventional dou- 
ble integral over coordinates in a manner that samples 
every pair (r, r') only once, rather than twice (with op- 
posite ordering, which requires the corrective factor of 
|). Each interaction term, defined by a pair (r, r') and 



one of the four combinations (c, a'), can then be decou- 
pled by a unique complex field 4> aa i (r, r'). But in order 
to write the result using the conventional double integral, 
we must create a copy of the Hubbard-Stratonovich field 
via </w (r, r') = -c?W(r',r): 



S sc = Jdtd\fl(D 

z<0 

+ 1 f dtd 3 rd 3 r'J2 



PI- 1 I 

2m 



(81) 



z,z'<0 



/\|2 



-|(^w(r, r')| 



V(r-r>) 
+ ^ a ,{v,v')fl{v)fl{v')+h.c 



Note that V(r — r') —> at large |r — r'| merely quenches 
the fluctuations of the corresponding (f) aa i (r, r') and pre- 
vents them from mediating long-range forces. This is 
also what keeps the action scale linearly with the system 
volume. The gauge transformations of 4> aa >(r,r') are: 



.(r,rO-> 



,(r,rV A(r) e lA(r,) 



(82) 



Next, we need to express the action in the Bogoliubov- 
de Gennes form in order to integrate out f a . This 
time, however, the presence of 4> aa >(r,r') makes the 
Bogoliubov-de Gennes Hamiltonian "HbcIG highly non- 
diagonal in the coordinate representation. Let us orga- 
nize all SC's fields into a vector F whose four-component 

blocks F(r) = f/ t (r), /|(r), / ; (r), /^(r)) are indexed by 

the coordinates r. Similarly, "HbcIG is a block-matrix 
whose blocks are indexed by coordinate pairs (r,r'): 



/ fi f D o - S " Msc) S(r' - r) 



6^Bd G (r,r') 



(r,r') 



b* n {v',v) 
80 



6 A) + ^+Msc M(r'-r) 



(r',r) 



80 



Z> n (r, r') 
80 



80 



r 



8 A-l-feU(r'-r) 



(r,r') 



(r',r) 



In order to represent the SC-2DES interface tunneling 
from ([67]) . we must also reorganize the 2DES's fields 
%jj Ta into a vector ^> compatible with F. *S> must be made 
of four-component blocks indexed by the coordinates r, 
and cannot have any leftover orbital dependence. These 
restrictions lead us to the blocks: 

*(r) =^y r (z) (ip Tt (x, y), ijil^x, y),ip Tl (x, y), ^(x, y) 

(84) 



is: 



S = J dtd 2 rip jf (D + H )iP 

2DES 



l - J did 3 rdV£ 



-(r,r')| 5 



2 

z : z'<0 
1 



F(r-r') 



2 F^H BdG F - - (^IF + h.c.) 



(85) 



The Bogoliubov-de Gennes form of the complete action The matrix X in the tunneling term is required by gauge 
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invariance, and its elements are given by (|68[) and 

^( r ) r ) = ^ & S(^ X x+yy,r''^xx+yy,r'^xx+yy,r'^xx+yy,r') ' 

Integrating out the SC's fields amounts to integrating 
out the Grassmann number components of the vector 
F. Since F groups together all available SC's Grass- 
mann fields, and not just a half of them as usual, this 
Gaussian integration produces the Pfaffian of the matrix 
%BdG; which we immediately relate to the determinant, 
Pf('HBdG) = \/det(%BdG)- The resulting effective theory 
of electrons inside the TI is: 



The generated quadratic coupling of VP in ([86)1 now takes 
the form $t (H s -H t )V, where $ = H' 1 !^. Since 
■H s contains no spinful triplets <j) aa , the components of W 
are operators that create or annihilate excitations with 
well-defined spin S z = ±|, despite being linear combi- 
nations of bare electron creation and annihilation opera- 
tors. Hence, we can think of 4> as the vector of quasipar- 
ticle field operators, which can further interact via the 
W (H s — Ht) \P term. Specifically, we are interested in 
their triplet-forming interactions: 



S e g = J dtd 2 r^ (D + H e )^ 



(86) 



2DES 



2 

z,z'<0 



dt d 3 r 



/■vl 2 



^W(r,r')| 
V(r-r') 



-i*ti^^ G it*_ itrlog^BdG ■ 

In order to learn something about the inter-orbital 
triplet pairing we need to analyze the inverse Bogoliubov- 
de Gennes matrix 7i^ G . To that end, let us separate 
Hb<ig = H s + Ht into pieces that couple electrons of op- 
posite and equal spins, H s and Ht respectively. Using 
the block-notation, 



H t (r,r') = 



80 
80 
80 

f>* n (r',r) 



H, 



HsdG 



4>n=4>U=0 



The H s part contains various singlet fields 4>aai where 
a = —a. Since <fi aS at two different locations (ri,r[) 
and (r2,r' 2 ) are not isolated by any physical symme- 
try, the effective action contains phase-locking couplings 
(jr (jS {Yi,r' l )(j) t7S (Y2,Y' 2 ) + h.c. that make all fields (j) GS si- 
multaneously condensed in the SC. The condensate am- 
plitudes are regulated by the microscopic interaction 
V(r — r') and decay rapidly when the separation between 
the end-points r and r' grows. At any rate, H s is char- 
acterized by a finite energy scale |^o| that corresponds 
to the superconducting order parameter. In contrast, 
the triplet fields 4> aa are not condensed in the SC. Their 
fluctuations at distances comparable to the quantum well 
thickness, |r — r'| sa d, have an effect on the triplet pair- 
ing in the TI, but the amplitude of these fluctuations can 
be considered small in comparison to \4>o\- Therefore, we 
may treat Ht as a perturbation to H s and write: 



^BdG — (H s +H t ) 1 — H s 1 ^ {-H t H s *) 



n=0 



H s 1 (H s —H t )H s l 



(88) 



|0 r l2 



z,z'<0 







(89) 



E 



^^^^(^(j'j + h. c . 



where r_L = (x, y) and we used the individual components 
ip Tcr (r±) of W defined by: 

This definition ensures that the quasiparticle operators 
ip Ta have the same dimensions, quantum numbers and 
normalization as the 2DES's electron operators ip Ta . The 
dynamics of uncondensed triplets <$> aa can be approxi- 
mated by a Gaussian: 



80 


80 1 


30 n (r,r') \ 


80 i 


W u (r',r) 


80 


u( r ' r ') 


80 


80 


80 


80 


80 J 



l - J d 3 rdVdt £ 

z,z'<0 



K CT (r,r')| 2 
V'(r - r'j 



(90) 



with a renormalized positive coupling V'. It is then 
(g^aight-forward to integrate out <p aa and obtain the in- 
duced interactions in the triplet channels: 



S% = -\ J dtd\d\'V'{v-v')Y, 



(91) 



*E 

{n} 



z,z'<0 



\<h\ 



x ^^rj^^)^^)^^) + • • • . 

Note that the effective action in the 2DES does not con- 
tain a direct order parameter coupling in the triplet chan- 
nel, because the SC has no triplet condensation. A direct 
coupling between 2DES's triplet pairs and SC's singlet 
pairs is forbidden by symmetry. 

We conclude that the Cooper pair tunneling mediates 
effective interactions among the TPs electrons that are 
short-ranged, attractive and weak if the SC's pairing gap 
|(^o| is large in comparison to the electron's SC-TI hop- 
ping energy scales t T (z sa 0). The induced pairing ten- 
dency at spatial separations r is related to the actual 
phonon-mediated attractive interaction V(r) that forms 
Cooper pairs in the SC. The induced interaction in the 
TI is generally capable of mixing the orbital content and 
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forming inter-orbital triplets as well as intra-orbital sin- 
glets. 

We must note that the strength of inter-orbital pair- 
ing in the quantum well 2DES, mediated by the virtual 
Cooper pair tunneling, depends on the renormalized SC's 
bulk interaction potential at electron separations compa- 
rable to the quantum well thickness, V'(d). This may be 
much smaller than the direct intra-orbital singlet pair- 
ing potential for two reasons. First of all, the quantum 
well may be a few crystal unit-cells thick, so that d is 
somewhat larger than the bulk SC's screening length l s , 
which is roughly equal to the inter-atomic separation. 
Even though d <~ l s by order of magnitude, the expo- 
nential suppression of the screened interaction strength 
with distance can make V'(d) weak. Second, the inter- 
action renormalization in the SC's bulk is contributed 
by the statistical (Pauli) repulsion between electrons of 
the same spin within the SC. This surely weakens V'(r) 
at separations r <~ d where electrons are quite delocal- 
ized. For these reasons, triplet pairing in the quantum 
well is most likely shaped by the direct SC's phonon 
mechanism, whose backbone Coulomb forces are spin- 
independent and felt across the quantum-well thickness 
without being jeopardized by screening, as we empha- 
sized earlier. In contrast, singlet pairing is contributed 
by the Cooper pair tunneling, but not strongly unless the 
hopping t T is unusually large or the SC order parameter 
|^o I is small. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we qualitatively analyzed the influence 
of a conventional superconductor on the electron dynam- 
ics in another non-superconducting material near their 
common interface. We identified two mechanisms of in- 
fluence, the coupling between the material's electrons and 
superconductor's phonons, and the virtual Cooper pair 
tunneling across the interface. We found that this yields 
induced effective interactions among the material's elec- 
trons near the interface, as well as a direct symmetry- 
breaking term in their effective action that constitutes 
the conventional proximity effect. 

The phonon proximity effect creates a retarded spin- 
independent attractive interaction that can in princi- 
ple be large enough to overcome the Coulomb repulsion 
and lead to independent Cooper pairing among the non- 
superconducting material's surface electrons. We calcu- 
lated the minimum electron-phonon coupling across the 
interface that leads to this proximity-induced pairing. Its 
value is related to the superconductor's electron-phonon 
coupling, the spatial extent d of the material's affected 
electrons near the surface, and the crystal's elasticity 
near the interface. Our generic conclusion is that any 
heterostructure with a sufficiently small d and thin in- 
terface will satisfy the pairing condition, by the virtue 
of the analogous condition being satisfied inside the bulk 
superconductor. 



If the superconductor's phonon-mediated pairing in- 
deed occurs among the material's surface electrons, then 
the critical temperature of the resulting two-dimensional 
superconductor is mostly independent of its surface den- 
sity. This does not hold at very low densities, near the 
surface superconductor-insulator transition. However, 
the character of the surface superconductivity changes. 
Cooper pairs can exist as two-body bound states in two 
dimensions, so their quantum phase transition is driven 
by quantum fluctuations, not pairing. These bound-state 
pairs are as large as the pairs in a dense two-dimensional 
BCS superconductor, so the effect of such fluctuations 
is visible only at extremely low densities where they do 
not overlap. More importantly, bound-state Cooper pairs 
can exist as gapped bosonic excitations in a (band) in- 
sulating state of the material's surface electrons, whose 
energy lies within the fermionic bandgap. Their binding 
energy is of the same order as the critical temperature of 
a dense BCS surface superconductor, and their gap can 
be controlled by a gate voltage. The existance and ac- 
cessibility of the "fluctuation regime" dominated by the 
dynamics of surface bound Cooper pairs is perhaps the 
most interesting perk of the interaction proximity effect 
because it can enable engineering novel correlated two- 
dimensional states of quantum matter. 

The tunneling mechanism, on the other hand, im- 
prints the superconductor's electron-electron interactions 
on the material's surface electrons via the superconduc- 
tor's Cooper pair dynamics. The induced interaction 
here is of the same kind as in the superconductor, i.e. 
short-ranged, attractive, and acting among electrons of 
any spin state. However, its strength is spin-dependent 
because the tunneling Cooper pairs are singlets. The 
induced surface interaction in the triplet channel is ex- 
pected to be weak, but it may nevertheless be help- 
ful to the phonon-mediated spin-independent interaction 
in stabilizing triplet Cooper pairs against the Coulomb 
repulsion. This is particularly relevant for materials 
that provide orbital electronic degrees of freedom, such 
as topological insulator quantum wells. They can host 
orbital-singlet spin-triplet Cooper pairs whose energy is 
either not prohibitively high, or even very competitive 
as a result of the spin-orbit coupling. Correlated quan- 
tum phases of triplet Cooper pairs in topological insula- 
tor quantum wells are a particularly attractive possibility 
because quantum fluctuations can in principle stabilize 
novel incompressible quantum liquids with non-Abelian 
quasiparticles. 

The Cooper pair tunneling is most effective in the sur- 
face singlet channels, which are directly coupled to the 
superconductor's order parameter. The effective action 
contains an explicit U(l) symmetry-breaking term as a 
result of this, the usually considered proximity effect. 
Such a term always produces singlet pairing correlations 
among the material's surface electrons, but their mag- 
nitude also depends on the bare surface electron's spec- 
trum. If the material's surface is normally metallic, the 
ensuing pairing correlations will be strong. This effect is 
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particularly interesting in the topologically protected he- 
lical metals on the surfaces of bulk topological insulators. 
The naively singlet pairing term actually gives rise to a 
very interesting surface superconducting state shaped by 
the spin-orbit coupling, whose topological defects bind 
Majorana quasiparticles. The analogous effect on band- 
insulating surfaces, i.e. topological insulator quantum 
wells, is diminished or absent. There, the pairing term is 
a small perturbation to the surface bandgap, so the ac- 
tual intrinsic singlet correlations are small. Instead, the 
main effect of Cooper pair tunneling are induced inter- 
actions in the singlet channels. They can shape bound- 
state Cooper pairs below the fermionic bandgap, whose 
dynamics can be made important by manipulating the 
surface chemical potential via the gate voltage. The 
ensuing singlet pairing is normally poised to win over 
triplets, unless the triplet dynamics becomes enhanced 
by the spin-orbit coupling. 

The main motivation for this and our previous 
wor k 38 ' 39 is the prospect of creating fractional incom- 
pressible quantum liquids in topological insulator quan- 
tum wells. The emerging picture is that there are no fun- 
damental obstacles to obtaining such novel strongly cor- 
related quantum states. In fact, the fundamental phys- 
ical principles that we explored theoretically encourage 
the existance of these states in the kind of systems we 
considered. The main question is if the microscopic pa- 
rameters of the available materials favor the existance 
of these states in devices that we can make, or conspire 
against them. This is a microscopic question that can be 
ultimately answered only by systematic experiments, or 
perhaps by ab-initio modeling of electronic spectra and 
their coupling to phonons in heterostructures. Neverthe- 
less, there are certain phenomenological conclusions that 
we could make so far. The spin-orbit coupling in topo- 
logical insulator quantum wells produces a fairly large 
SU(2) cyclotron energy scaled 8 ., which can be of the or- 
der of 100 meV. A crude calculation in this paper shows 
that the proximity-induced phonon-mediated interaction 
in the quantum well alone could be large enough for the 
two-dimensional pairing to occur in certain heterostruc- 
tures. These quantum wells are naturally insulating, so 
if they can be tuned through a phase transition by a gate 
voltage, their density of mobile spinful particles, triplet 
Cooper pairs, can be made small and comparable to the 
SU(2) spin-orbit flux-quantum density 3 - 8 -. The quantum 
fluctuations associated with Cooper pairs are large in this 
regime, so most essential ingredients for fractional incom- 
pressible liquids are present. 

What we cannot predict is whether the quantum wells 
can be made clean enough to obtain sufficiently high mo- 
bility (important for fractionalization) and spatially flat 
potential landscape (important for the gate voltage tun- 
ability) . We also cannot predict exactly how large the net 
proximity-induced pairing interactions are in the quan- 
tum well. They have to be large enough in order to ensure 
pairing. The required electron-phonon coupling across 
the interface can be maximized by eliminating the tun- 



neling barrier in the FigQ] and keeping the quantum well 
thickness d very small. However, this simultaneously di- 
minishes ones ability to control the chemical potential in 
the quantum well via the gate voltage. It might be neces- 
sary to find an optimal barrier and quantum well thick- 
ness that achieve a compromise between the electron- 
phonon coupling and tunability. The choice of materials 
might also be important, since their interface work func- 
tions can make it harder or easier to place the chemical 
potential in the quantum well near its conduction (or va- 
lence) band where bound-state Cooper pair modes should 
exist. 

There are many other kinds of heterostructures and 
interfaces that we did not address in this paper. One 
obvious question is whether the proximity effect created 
by a high-temperature superconductor could give corre- 
lated states of matter a better chance. Certainly the con- 
ventional proximity effect caused by the U(l) symmetry 
breaking would be stronger. However, the induced inter- 
actions need not be stronger at all. We have found, in 
fact, a certain inverse relationship between the supercon- 
ductor's order parameter magnitude (critical tempera- 
ture) and the strength of proximity-induced interactions, 
because it is the quantum fluctuations that generate the 
effective interactions. Furthermore, the cuprate super- 
conductors have a pairing glue whose not-yet-understood 
dynamics is visible only at very short length-scales and 
hardly transferable to the other material's surface elec- 
trons. In the simple picture of spin-exchange-mediated 
pairing in cuprates, there is no mechanism for the mate- 
rial's surface electrons to experience the superconductor's 
spin exchange across the interface (the material's surface 
is not close to being a half- field fermionic Mott insulator). 
This eliminates the possibility of the proximity-induced 
triplet pairing, additionally aggravated by the strongly 
spin-dependent nature of interactions in cuprates, which 
is passed on to the material's electrons by tunneling. 

Other important proximity effects are created by mag- 
netic rather than superconducting materials. Here, too, 
we find no significant analogue of the phonon mecha- 
nism. Since magnetic orders are stabilized by Coulomb 
interactions, photons would be the proper substitute for 
phonons in magnetic proximity effects. The ensuing mag- 
netic interactions are dipole-dipole, typically not very im- 
portant at short length-scales. The effective surface in- 
teractions mediated by electron tunneling can be deduced 
using the similar approach as in the second part of this 
paper, but done in the particle-hole channels. We expect 
a similar phenomenological outcome, that the magnetic 
material's repulsive electron-electron interactions are im- 
printed in a weaker form on the surface electrons across 
the interface. There is, of course, an induced Zeeman 
coupling as well, the result of symmetry breaking usually 
considered as the proximity effect. The induced interac- 
tions could in principle lead to new surface phases, es- 
pecially if the surface dynamics creates opportunities for 
instabilities in the presence of interactions that do not 
exist in the magnetic material. 
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better one) to label a time variable in this appendix; this 
should not be confused with the orbital labels from the 
rest of the paper. 

The two-body Green's function is: 

x^(^-f)^(^ + ^)|0), (Al) 



Appendix A: Bethe-Salpeter equation 

Here we review the properties of the two-particle 
Green's function Q that appears in the Bethe-Salpeter 
equation (f29|) . Our specific goal is to establish the rela- 
tionships (HU), ([25]) and (1321) that led us to the formu- 
lation of the Schrodingcr equation for the Cooper pair 
bound-states. Note that we use the symbol r (lacking a 



where ,ip are the single-particle creation and annihi- 
lation operators in the Heisenberg picture, and X^ — 
(R, r) are the center-of-mass and x^ — (r, t) the rela- 
tive coordinates of the two particles. After expanding 
the time-ordering, partially applying the time evolution 
■0(r ± t/2) = e lHT '4){±t/2)e- tHT and inserting the iden- 
tity operator resolved in terms of all exact eigenstates \n) 
of the many-body Hamiltonian, we obtain: 



E p e-^l r - r V p (R,r,t)J-t(R',r',t') , r - H > r' + ^ 
T, h e- iEhlr - T 'lF h (R,r,t)7l(R',r',i') , r + M < r' - M 



(A2) 



where we emphasized only the contributions of two- 
particle \p) and two-hole \h) excitations, and defined 

J" p (R,r,t) = (0|TV ( R + V ' ( R ~ lp) 

J- fe (R, r ,i) = <W> ( R + 1^)^(^-1,-0 |0> • 

E p and Eh are the exact energies of the two-particle and 
two-hole eigenstates. Taking the Fourier transform over 
r yields: 

G(R, ft, r, t; R', ft', r', t') = (2n)8(tt - ft') x (A3) 

E 

L P 

E 



Jp (R, r, t)Ti(R', r', t') i(n _ EB+i0+) m+l^j _ 
ft - £p + i0+ 

■F fe (R, r, t) j|(R', r', t') i(a+Eh _ i0+) m+£± 

n + E h -io+ 



We will be interested in a particular two-particle bound- 
state \p) with energy E p — £ , and set accordingly ft 
arbitrarily close (or at) £: 

g( R,n,r,,;R/,n', r -,0 ^ (R ^"^ r ' i,) 

ft — e + z0 + 

x (27r)<5(ft - ft') . (A4) 

The correction to this expression is what we called the in- 
coherent part in (|24l) . Note that the explicit exponential 



factor involving \t \ + \ t'\ is absent in the written coherent 
part. We may next carry out the Fourier transform with 
respect to t, t 1 : 

ft — t + iQ + 
x (27r)<5(ft - ft') , (A5) 



where 



J"(R,r,w) = j dh 



(A6) 



Note that 

/j 
— ^(R,r,«)= J-(R,r,i = 0) (A7) 

= (O|^(R+0V(R-|) |p)s*(R,r) 

is the standard stationary two-body wavefunction of the 
bound-state with the center-of-mass at R and the rela- 
tive displacement between particles r. This leads to the 
equation (|25)l . In the usual circumstances, we can sepa- 
rate variables ^(R, r) = 5' cm (R)-0(r), where ip(r) is the 
wavefunction of the relative motion of the two particles. 
Lastly, we may carry out the Fourier transform to the 
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momentum space: 
where 

J-(P, p,u) = J d 2 R d 2 r e- l(PR+pr) J-(R, r, u). (A8) 

Note that the total energy £ depends on the total con- 
served momentum P conjugate to the center-of-mass po- 
sition R. 

The Bethe-Salpeter equation (|29p is obtained by sub- 
stituting (| A8|) into the Dyson equation (|23l) . taking the 
f2 — > £ limit and neglecting the incoherent parts of the 



Green's function. Specifically, the bare Green's function 
Go that appears as a separate additive term in the Dyson 
equation is negligible at the bound-state pole of the exact 
Green's function Q. After this term is neglected, the com- 
mon factors of from the exact Green's function can be 
canceled out since nothing operates on them. Finally, by 
the separation of variables .F(P,p,w) = J-" cm (P)F(p, to) 
into the center-of-mass and relative coordinates, we can 
cancel out the common factors of J-" cm (P) to obtain (|29l) . 

The function J^P,^) is still a Green's function with a 
frequency dependence. We can try to elucidate the mean- 
ing of this frequency by the same approach that we took 
to reveal the meaning of the center-of-mass frequency 
ft. To that end, we will carry out the time-ordering a bit 
more carefully than above (assuming fermionic particles): 



JYP,^)=E 



JHt/2 



V(R+§). 
R 



iHt/2^ 

1*1 



e > Ht / 2 \m)(m\e- lHt / 2 ^(R- 
5j e iHt / 2 \m)(m\e iHt / 2 ip (R + 

(0|V(R+§) |m)H^(R-|) \p) 
-<0|V(R-§) |m)(m|V(R+§) \p) 



iHt l 2 \p) , t>0 



t > 
t < 



(A9) 



The sums run over all exact single-particle excitations 
|m). The Fourier transform yields: 



(o|V(R- 



i)(m\ip (R- 



\p) 



+ 



-25, 

(0|^(R-§) |m)(m|V(R+§) \p) 



i0+ 



-E m \- iO 



(A10) 



If the two-body spectrum contains a bound state \p), 
then all single-particle states |m) have energies E m that 
lie above the bound-state energy E p by a finite amount, 



E p — 2E m < 0. The largest E p — 2E m is the binding en- 
ergy e < 0, and setting w = ±e/2 should generate a pole. 
In the absence of a bound-state, the would-be pole is em- 
bedded in the two-particle continuum and smeared out 
by the integral (sum) over |m). In a conventional bulk 
BCS superconductor, however, the pole is resurrected de- 
spite the absence of a Cooper pair bound-state due to the 
pairing gap A : E p = 0,E m = y/E 2 (p) + Ag, implying 
E p /2 — E m < for all |m). Clearly, the pairing gap and 
a Cooper pair binding energy (when applicable) play the 
same physical role. The function F(p,u) turns into the 
anomalous propagator in the BCS theory. 
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